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ABSTRACT 


This study utilizes che RFLAP4&/MODS Computer 
Code to analyze the effects of critical flow 
modeling on the thermal-hydraulte transient 
response of a General Flectric bhbottiinp water 
reactor to a major primary coolant line rupture. 
Included itn the study 13 a presentation of the 
equations, assumptions, and limitations of the 
critical flow models availahle for use in RFIAP4. 
Additionally, an evaluatton of a temporary 
solution to a RELAP& coding error associated with 
stagnation properties calculation {is presented. 
The results of this investigation indicate 
that; (1) A solution to the stagnation properties 
calculational error that extats in the RFLAPG 
ode, when applied to the Fvaluation Mosel, 
(PELAP4-FM) provides no conservative evaluation, 
relative to the Standard Model (RFLAP4&-SM), 
regacdlesan of which of the five critical flow 
models avaitable in the code in selected; (2) Ge 
the five critical flow models avatiahble for use itn 
RRLAP4, the Moody model and the HWenry-Fauake model 


are nearly equivalent in their relative cepreen of 
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Mementua/Htomogeneounr Foutlibriun model, the 


Hose pencout Fawilibedum model, and the fonic 


models; (3) In order to best alleviate the 
stagnation properties coding error, the flow areas 
for the volumes immediately upstream from the 
recirculation line rupture should the increased 
such that the area ratio (break flow area / 


upstream volume flow area) is 0.7. 
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ABSTRACT 


This study utilizes the RELAP4/MODS5 Computer 
Code to analyze the effects of critical flow 
modeling on the thermal-hydraulic transient 
response of a General Electric boiling water 
reactor to a major primary coolant line rupture. 
Included in the study is a presentation of the 
equations, assumptions, and limitations of the 
critical flow models available for use in RELAP4. 
Additionally, an evaluation of a temporary 
solution to a RELAP4 coding error associated with 
stagnation properties calculation is presented. 

The results of this investigation indicate 
that; (1) A solution to the stagnation properties 
calculational error that exists in the RELAP4 
code, when applied to the Evaluation Model, 
(RELAP4-EM) provides a conservative evaluation, 
relative to the Standard Model (RELAP4-SM), 
regardless of which of the five critical flow 
models available in the code is selected; (2) Of 
the five critical flow models available for use in 
RELAP4, the Moody model and the Henry-Fauske model 
are nearly equivalent in their relative degrees of 
conservatisme Following these models, in order of 
relative degree of conservatism, are the Modified 
Momentum/Homogeneous Equilibrium model, the 


Homogeneous Equilibrium model, and the Sonic 





model; (3) In order to best alleviate the 
stagnation properties coding error, the flow areas 
for the volumes immediately upstream from the 
recirculation line rupture should be increased 
such that the area ratio (break flow area / 


upstream volume flow area) is 0.7. 
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void fraction 
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time step size 
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mass flow rate per unit area (mass flux) 
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enthalpy 
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coefficient of isothermal compressibility 
lumped parameter accounting for friction 
and momentum changes 
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entropy 








surface tension 
velocity 
specific volume 
mass flow rate 
quality 


distance in direction of fluid travel 


SUBSCRIPT NOTATION 


t+ At 


UP 


downstream or back 

critical 

liquid 

vapor 

junction 

Saturated conditions 
transition completion point 
time at particular time step 
time at following time step 
upstream 


Stagnation 
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Ll. INTRODUCTION 


1.1 GENERAL 


The United States, in its fourth decade of nuclear 
reactor research and operation, now has seventy-two civilian 
power reactors in operation and has’ generated over a 
trillion megawatt-hours of electricity. Since the inception 
of the nuclear power industry, safety has been the major 
consideration in the design, construction, and operation of 
nuclear power plants. 

Modern experimental programs in water-reactor safety 
are directed toward improving the ability to understand 
specific phenomena and the results of tests, with the goal 
of applying this understanding to the interpretation of 
"similar" phenomena that might take place if ae reactor 
accident were to occure Complex digital computer codes have 
been developed to predict the thermal-hydraulic transient 
response of nuclear reactor systems to hypothesized 
accidents. Verification of specific portions of these 
computer codes is the goal of an ongoing experimental water- 


reactor safety program. 





1.2 LOSS-OF-COOLANT ACCIDENT (LOCA) 


A LOCA is) generally defined as any accident which 
causes the coolant to be lost from the cooling system 
serving the reactor core. The design-basis LOCA is. that 
primary system pipe rupture which results in maximizing the 
fuel-rod cladding temperatures attained during the accident. 

For boiling water reactors (BWRs) the design-hasis LOCA 
is a double-ended guillotine break of one of the two coolant 
recirculation lines. The coolant recirculation lines are 
the primary means of maintaining coolant flow through the 
reactor core, and are the largest diameter piping in the 
primary system. Emergency core cooling systems (ECCS) are 
designed to provide cooling to the core in order to prevent 
fuel element melt down in the event of a design-basis LOCA. 
Future references to LOCA will imply design-basis LOCA 
unless otherwise indicated. 

Within the present regulatory framework, the thermal- 
hydraulics of a LOCA can be described as’ having’ three 
distinct phases known as blowdown, refill, and reflood (2). 
(1) Reactor Blowdown: During the reactor depressurization or 
blowdown phase, reduction of core flow and increase in the 
local quality result in departure from nucleate boiling 
(DNB) and dry-out of the core- Following the establishment 
of the conditions causing DNB, the heat transferred from the 


fuel rods to the coolant decreases markedly, and_ the 





temperature of the fuel rods, and the fuel-rod cladding, 
rises rapidly. Depressurization continues until primary 
system and containment pressures are equalized and the break 
flow is essentially zero. 
(2) Refill and Reflood: At the end of the blowdown phase, 
the liquid coolant inventory in the system is small, and the 
core may be completely immersed in steam. The system 
pressure is low, on the order of 50 to 100 PSIA. As the 
ECCS starts injecting emergency coolant, refilling of the 
lower plenum begins- For conservative safety evaluations, 
heat transfer to the coolant is often totally tgnored during 
this refill phase- When the liquid level reaches the hottom 
of the core, the reflooding phase beginse- As_ the liquid 
level rises in the core a tremendous amount of heat is 
removed from the fuel rods due to vaporization of the 
incoming liquid. Peak cladding temperatures will he 
attained during this phase of the accident. 

For an excellent, detailed physical description of a 
BWR LOCA the interested reader is referred to references (3) 


and (4). 


1.3 CODE DESCRIPTION 


RELAP4 is a computer program, written in Fortran IV, 


that was developed primarily to describe the thermal- 


hydraulic transient behavior of water-cooled nuclear 





reactors subjected to postulated accidents such as’~ those 
resulting from loss of coolant, pump failure or nuclear 
power excursions. Fundamental assumptions inherent in the 
thermal-hydraulic equations used in the code are that a two- 
phase fluid is homogeneous and that the phases are _ in 
thermal equilibrium. Models are available to modify the 
homogeneous asSumption. 

The program is currently available in its fourth major 
revision, hence RELAP4. The edition of the program 
installed at The Pennsylvania State University is 
RELAP4/MODS5S-UPDATE2 (henceforth referred to as RELAP4G). 

There are four major options available in RELAP4. The 
Evaluation Model option, hereafter referred to as RELAP4<-EM, 
is intended to provide conservative results which meet most 
of the requirements for licencing under Appendix K to 10 CFR 
Part 50. Another option available, called the Standard 
Model (RELAP4-SM), is used to provide the best estimate of 
the actual expected response during a LOCA. While RELAP4 
was primarily designed to be a blowdown code, it is capable 
of representing the LOCA transient through reflood in 
pressurized water reactors (RELAP4=-FLOOD option). Also 
available is the containment option (RELAP4-CONTAINMENT), 
which can be used to evaluate the effectiveness of ice- 
condenser containment systems- The two main options that 


will be used in this study are RELAP4-EM and RELAP4~SM. 





RELAP4 requires numerical input data which completely 
describes the initial conditions and the geometry of the 
system being analyzed. Included in the input data _ are 
physical characteristics such as fluid volume, geometry, 
pump characteristics, power generation, and material 
composition. Transients can be initiated by the control 
action inputs to the programe These controls can describe 
breaks in fluid piping, core power level variations 
including a SCRAM, valve actions, etc. For each time 
advancement the program computes fluid conditions’ such as 
flow, pressure, quality, and mass inventory. Also computed 
are thermal conditions within the solid materials’ such as 
temperature profiles and power generated, and the fluid- 
solid interface conditions such as heat flux and surface 
temperature. 

A complete theoretical development of the code, as Well 
as implementation procedures, i185 contained in the 
RELAP4/MOD5 Users Manual (5)- An excellent description with 
specific reference to BWR LOCA analysis is also available in 
a masters thesis by Laird (3). 

The reactor system used for this study is the Hope 
Creek Nuclear Power Generating Station. This General 
Electric BWR was modeled for RELAP4 by Idaho National 
Engineering Laboratory (INEL) personnel, and was included as 


a sample problem with the RELAP4 transmittal tape. 





The thermal-hydraulic model of Hope Creek, depicted in 
Figure 1, consists of twenty-three control volumes, thirty- 
three junctions, and twenty-five heat slabs. Table 1, pages 
8 and 9, presents a description of the twenty-three control 
volumes used in this’ model. Junctions simply connect 
control volumes, and heat slabs represent sections’ of the 
physical system which can absorb, radiate, store, or 
generate heat. RELAP4-EM and RELAP4-SM data input decks 


used in this study are discussed in Appendix A. 


1.4 CRITICAL FLOW 


"When a fluid expands from a compressed state to 
arbitrary ambient (or "receiver") conditions tin passage 
through an outlet, the flow rate is always less than a 
certain maximum or "critical" value. After this critical 
condition has been reached, further reductions in receiver 
pressure leave the flow rate unaltered, serving only to form 
a steep pressure gradient at some location in the outlet 
passage (6)." This condition is commonly called critical, 
or choked, flow. 

At the instant of recirculation line rupture the flow 
rate out of the break increases until critical flow 
conditions are reached. There are five critical flow models 
available in RELAP4 to predict this critical flow rate. 


They are, (1) The Sonic model, (2) The Moody model, (3) The 
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Table 


No. 


lL. Volume Descriptions in RELAP4 Mode] 
Description 
Upper Plenum and Steam Separators 


Upper Downcomer (This volume extends from the feedwater 
entry point to the bottom of steam dryers.) 


Steam Dome (This volume contains the high quality steam 
located above the dryers.) 


Lower Downcomer (This volume contains the subcooled water 
which enters the downcomer from the feedwater entry 
point. The lower downcomer describes the region of the 
downcomer between the feedwater entry point and the 

entry point of the feedwater into the lower plenum.) 


Broken Loop Recirculation Suction Line (This volume 
models the 26" pipe which pulls the saturated water 
from the lower downcomer into the jet pump of the 
broken loop. Although the break occurs along this suc- 
tion line, the actual break volume is modeled by 

volume 21.) 


Broken Loop Recirculation Pump 

Broken Loop Recirculation Discharge Line (This volume 
models the 24" line which returns the water from the 

recirculation pump to the jet pump intake.) 

Intact Loop Recirculation Suction Line 

Intact Loop Recirculation Pump 

Intact Loop Recirculation Discharge Line 

Lower Plenum 

Heated Core Sections (With the addition of volume 23, 
these volumes model the core. To achieve relatively 

equal thermodynamic characteristics, the heated core 

sections vary in height from 2.85 feet to 1.00 feet.) 
Guide Tubes (The control rod guide tubes in the lower 


plenum and the unheated segment of the core are modeled 
as one control volume. ) 





Volume No. 


18 


19 


20 


Zi 


22 


235 


Table 1. (cont.) 
Description 


Core Bypass (This volume represents the region between 
the fuel assemblies in the core.) 


Intact Loop Jet Pump (This volume models the intact loop 
jet pumps from the intake to the exhaust into the lower 
downcomer. Each bank of ten jet pumps is modeled as 

one volume. ) 


Broken Loop Jet Pump (This single volume models the 
bank of ten jet pumps in the broken loop.) 


Break Volume (This volume is specified between the 
broken loop recirculation discharge line and the lower 
downcomer. ) 


Containment 


Unheated Segment of the Core 
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Henry~Fauske model, (4) The Homogeneous Equilibrium model, 
and, (5) The Modified Momentum/Homogeneous Equilibrium 
model. These models are discussed in detail in Chapter Two 
of this study. 

The simple momentum equation for inertial flow rate is 
quite accurate at low flow rates but becomes”) greatly 
exaggerated as the driving pressure differential increases 
and critical conditions are approached. RELAP4 computes the 
flow rate using the simple momentum equation and compares 
this flow rate to that predicted by the’ selected critical 
flow model-~ The minimum is then selected as depicted by the 


solid line in Figure 2 (7). 


Critical Flow Criteria 
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FIGURE 2 
JUNCTION MASS FLOW RATE SELECTION 
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1.5 STAGNATION PROPERTIES 


All of the critical flow models available in RELAP4, 
with the exception of the Sonic model, require _ that 
stagnation properties be calculated for the volume upstream 
of the junction undergoing critical flow. These 
calculations are performed by a subroutine called STAGP. 

The stagnation pressure and enthalpy are calculated 
from the static pressure and enthalpy using a calculated 
Fluid velocity. An isentropic path from the static to the 
stagnation enthalpy=<-pressure point is followed. 

The stagnation enthalpy is calculated from the kinetic 


energy relationship: 


2 
hl wh, + Vo/2g 5 (1) 


where V is calculated from the following relationship: 
VeWv/A . (2) 
The quantities hy» W, and v are evaluated at the junction, A 
is the upstream volume flow area. 
In equation (1), V is not allowed to be greater than 
the isentropic sonic velocity of the upstream volume. 
Applying the first law of thermodynamics in 


differential form for a fluid undergoing a flow change 


provides, (Gibbs equation); 





bz 


Tds = dh - vdP . (3) 
For an isentropic process, ds = 0. Then, 
dP = dh/v. (4) 


Integrating this for the change in enthalpy and 


applying static and stagnation limits provides: 


h 
ee 2: a dh/v (5) 
1 
where v, the specific volume, is a function of P and h, and 
aa and Py are the stagnation and static pressures, 
respectively. 

Equation (5) is advanced from (hi »P,) to (hosP) by a 
fourth-order Runge-Kutta subroutine. 

It had been reported that there existed a logic problem 
in the stagnation properties subroutine, STAGP, requiring 
the elimination of this routine from the RELAP4 logic (8, 
9). It has since been discovered, and disclosed at the INEL 
Code Users Workshop conducted at INEL, Idaho Falls, Idaho 
8/28/78 thru 9/1/78, that the stagnation properties 
calculation problem was not due to a logic error in the 
STAGP subroutine, but rather was due to the fact that the 
code was predicting velocities greater than the isentropic 
sonic velocity for the volume immediately upstream of the 


break junction. 








i3 


As pointed out earlier, the code does not allow the 
calculated velocity in equation (1) to exceed the isentropic 
sonic velocity of the upstream volume. If the junction 
conditions are such that equation (2) predicts a velocity 
greater than sonic, erroneous values for the stagnation 
pressure and enthalpy will result. 

The "fix" suggested at the workshop was to artificially 
increase the upstream volume flow area (the denominator of 
equation (2)) in order to prevent’ supersonic velocity 
predictions from occurringe This flow area increase should 
be carried out for all volumes where critical flow is likely 


to occuLrs 


1.6 PURPOSE 


The purpose of this study is threefold; 

(1) Review the derivation of the equations for each of 
the critical flow model options available in RELAP4. 

(2) Determine by what amount the upstream volume flow 
area must be increased in order to alleviate the problem in 
calculating stagnation properties that was reported in the 
previous section. 

(3) Present a sensitivity study on the critical flow 
models available in RELAP4 and rate the models as to their 


relative degrees of conservatism. 
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2. CRITICAL FLOW MODELS 


2-1 INTRODUCTION 


As previously reported, there are five critical flow 
models available in RELAP4. They are: 

(1) The Sonic model (10), 

(2) The Moody model (11), 

(3) The Henry-Fauske model (12), 

(4) The Homogeneous Equilibrium model (13), and 

(5) The Modified Momentum/Homogeneous Equilibrium model 

(14). 

For each junction in the system the user must specify 
whether the code is to use one of these models or the simple 
momentum (inertial flow) equations to estimate the flow rate 
through that junction- In addition the user may select from 
eighteen possible combinations of these models.~ For example 
the user may elect to use the Sonic model when the fluid is 
subcooled and the Moody model when the fluid is saturated. 
These combinations, as well as the selection methods, are 
discussed in the RELAP4/MOD5 Users Manual (15). 

The purpose of this chapter is to present each of the 
critical flow models listed above, providing the basic 
assumptions and limitations of each. In addition, at the 
end of this chapter a brief description of the logic 


involved in junction flow rate estimation is given. 
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2-2 THE SONIC MODEL (10) 


As the name implies, this model assumes’ that’ the 
junction flow rate is that which will result in sonic 
velocity at the junction. 

The sonic velocity of the fluid is calculated fron, 


2 2 Fixe ’ , 
a= eae) gt (P ) /{X{C,, J+TP Met Kaeo ee 


g 
(6) 


+ (1-X)[C,, Se ees Vv, OP K 


= eee) ly 
£ f 


f 
where P’ signifies the derivative of the partial pressure 
with respect to temperature. 

It is important to note that the independent variables 
in equation (6) are based upon the upstream volume 
conditions. That is, the steam tables are entered with 
upstream volume specific internal energy and specific volume 
as the independent variables- All other upstream volume 
parameters, including pressure and sonic velocity, are 
dependent on these two quantities. The result is that the 
upstream volume critical velocity is calculated and then 
impressed upon the local downstream junction. 

The critical mass flow rate through the junction is 


calculated from the one-dimensional continuity equation, 


W= OA a ° (7) 
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The density used in equation (7) is ecalculated from the 
upstream volume density as follows. The upstream volume 
density is first adjusted for the frictional loss_ and 
kinetic energy changes to the junction. The fluid is then 
expanded isentropically from the Mach number just upstream 
from the junction to that in the junction, which is Mach 
one, when sonically choked. It should be noted _ that a is 
the flow area at the junction. 

The assumptions implicit in the use of this model are: 

(1) The critical flow rate is assumed to be that which 
will result in sonic velocity at the junction, 

(2) The critical flow rate is given by the one- 
dimensional continuity equation, 

(3) The junction critical velocity is based upon 
upstream volume fluid conditions, and 

(4) In calculating the density used in equation (7) 
isentropic expansion is assumed. 

This model is also capable of predicting critical flow 
rates when the upstream volume contains air or a mixture of 
air and saturated or superheated steam. This capability is 
relevant only to the containment option (RELAP4-CONTAINMENT ) 


and will not be further discussed. 





ey, 


2.3 THE MOODY MODEL (11) 


The Moody model is applicable only in the saturated 
region, and therefore is used primarily in low pressure 
systemse It is the only critical flow model available in 
RELAP4 which allows for slip between phases. In addition it 
does not rely on empirical data for its derivation. 

The derivation of the equation for critical mass flux 
starts with the continuity and energy equations for two- 


phase, annular flow; 
G=W/A = aV /Xv #® (1 -aQ@)V_/(1 = X)v (&) 
gx gZ f f 
rp EE Os CI ay ey a Un ere ere ay Be (9) 
oO g g, c f f c 
X and a@ are defined by, 
X = egal : l1=- X «= We/W (10) 
= e - Q es 
O AL/A : l A./A (11) 


The ratio of average vapor velocity to average liquid 


velocity is the slip ratio, §S ; 


S = Ue : (12) 
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Rewriting equation (8) and substituting equation (12) 


yields, 


UM = (1 - ove cael - X)ve = S e (13) 


Solving for the void fraction in terms of the slip 


ratio yields, 
a= 1/{1 + [Svp(1 - X)/Xv 1) (14) 


Assuming that the expansion is isentropic (s = s) and 


from the definition of quality, 
X sx = 
(s | S,)/8., (15) 
or 
l1=-X = (s, = Solu Ger ° (16) 
An expression for the mass flux G is obtained by 
combining equations (8) through (16) as follows; 


Rewriting equation (9) yields, 


2g Jth > - he = X(ho = he] = Welpe < ve (ix) /V2} Gi) 
g g 





WS) 


Substituting equation (12) into equation (17) and 
hat h =h -h lelds 
noting tha fs 8 f yie y 


2 2 
22 J Ch - he - re = Vale + (1 - X)/S°J ; (18) 


Rewriting equation (14) as, 
Q = aOR 4 Sv, (1 - X)]) ; (19) 


solving equation (8) for Ve and substitutinp equation (19) 


yields, 
Ve=2# G[Xv + Sv,(l - X)]J ° (20) 
g g f 


Substituting equation (20) into equation (18) = and 


solving for G yields, 


2 
G 2a dh, - he - Xhe,)/ 
CZ") 


({Xv, + Svp(1 - x))7tx + (1 - x)/8*}} 


Substituting equations (15) and (16) into equation (21) 


yields, 





2 


2 


G° ize Jth, - h, - bh, (es. - Bed eee! 


f fg 


2 
eae - arg a + vpecere - Spar C225 


alban: 


[(s) - ee + SS = s)/S’s 

Note from equation (22) that if h and So are known, 
then G is a function of S and P (ie. knowing the pressure 
and the fact that saturation exists fixes all of the 
quantities in equation (22) except ho» So? and S). 

Therefore, if G has amaximum it must satisfy the 


necessary conditions, 


a 
Q 


| 


(23) 


Q? 
Y 
rd 


and, 


rob 
O 


| 
tt 
Oo 


(24) 


a 
a=) 
YI 


assuming that S and P are independent. 
Applying equation (23) to equation (22) and solving for 
S yields, 


stBv fv y!/3 . (25) 
g f 


Thus for maximum flow rate S depends only on P. 
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When equations (25) and (22) are combined, a maximum G 
must satisfy equation (24)- Equations (22), (25), and the 


state equation, 


h «® h +h (s - s. )/s (26) 


were programmed for saturated steam/water properties (16) 
and solved by digital computer for values of pressure which 
satisfied equation (24) for input values of P and ih =. 
Calculations showed that G has a single maximum value for 
known Ee and ho? The maximum mass flux was’ thereby 
determined. 

The results of the above computations are available in 
tabular form for use in RELAP4. The tables are entered with 
upstream volume stagnation pressure and junction stagnation 
enthalpy and a single mass flux is returned. 

The Moody model is based upon an annular, two-phase, 
one~dimensional model with uniform axial velocity of each 
phase and thermodynamic equilibrium between phases. 
Additional assumptions are: 

(1) Both phases experience the same local. static 
pressure, 

(2) The flow is isentropic from entrance to exit, hence 


Stagnation enthalpy is constant, 


(3) The liquid phase is incompressible, and 
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(4) The slip ratio and pressure are independent 
variables. 

The Moody model is strictly applicable only in the 
saturated regime with stagnation pressure limits of DPS PA 
to 3000 PSIA. If the Moody model is used in the suhcooled 
regime the code will enter the critical flow tables at zero 


quality. 


2-4 THE HENRY= FAUSKE MODEL (12) 


The Henry~Fauske model was developed primarily to 
account for the nonequilibrium nature of fluid flow under 
critical conditions and is applicable for fluids in the 
subcooled as well as the saturated and superheated states. 

The derivation starts with the steady-state, one- 
dimensional continuity and momentum equations for two-phase 
flow; 


Liquid continuity, 


W.v, = A.V (27) 


Wv 8 AV (28) 





Zs 


Momentun, 


AdP sh CE + W, p> /8. dF, 23) 

For high velocity flows the wall shear forces are 

negligible compared to the momentum and pressure pradient 

terms. Therefore, with the substitution of the liquid and 
f 


vapor flow rate definitions, fue = X W and W = (1 - X)W), 


equation (29) becomes, 


dP Gd [X . (1 ) p)/8. (30) 
It is assumed, for fixed upstream Stagnation 
conditions, that X, ve Vee oe and v are either constant 


or composite functions of P and ze Hence, equation (30) 


becomes, 


] 


- d 
Pitan’, 105 Ane Seth FEELS: (31) 


j j 


where the subscript j indicates that the enclosed quantities 
are evaluated at the junction- At critical flow, the mass 
flow rate exhibits a maximum with respect to junction 


pressure; 
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Equation (32) can be applied to equation (31) to give 


an expression for the critical flow rate; 


-2 


-G., 


d 
= {ap (Xs + (1-X)IEC1-X)Sv_ + AP er (33) 


where, as before, S is defined by 5S = Ue Nes By taking the 


implied derivative, equation (33) can he expanded to 


-677 = (X[14X(S-1)]dv /aP + S[1 + K($-2)-X°(S-1)]dv,/dP 


+ vee + 2X(S-1l) + 2Sv,(X-1) + S$(1-2X)JdX/dP (34) 


+ X(1-X)[v,8 - v,/S1dS/dP} /g 8 - 


j 


Equation (34) can be simplified by using the following 
assumptions and relationships: 

(1) The vapor and liquid velocities are equal (no 
slip); 

S = 1 (35) 


(2) The liquid phase is incompressible; 


dv /dP = 0 (36) 





ILD) 


(3) The vapor behavior can be described by a polytropic 


process such that, 
dv /dP = v /nP (37) 
g g 


where n is the thermal equilibrium polytropic exponent 


derived by Tangren et al. (17) and given by, 


n= {(1 - X)Cy ic, + XJ/{((1 - X)Cy /C 


+ X/y ] (37a) 
f f : 


V 


(4) Negligible interphase mass transfer; 

X = X (38) 
(5) The empirical result; 

qs/dP], =a 0 . (39) 


Equations (35) through (39) serve to reduce equation 


(34) to 


-1/2_ 1/2 
Bo 


Ge els = — aN )dX/dP), 


O 


(40) 


Formulating the actual junction fluid quality in terms 
of the equilibrium quality at that location and making use 
of empirical data as well as some mathematical manipulation, 


equation (40) becomes 





2 6 


1/2 
= P + 
G | 8. { Xv, /0 


fe sds, /dP - (41) 


(v - v. ){{ (1 - X_)N/s 
B Oo w E E 


f 


-1/2 


{XC /Ps j 


P 


f POl in. = 177) 15) 
& 
g O 


The nonequilibrium parameter N is of primary import in 
the low quality regime. Based on experimental data (18), 
this parameter can be related to the equilibrium quality 


(X_ ) as follows, 


Py 


Xf Odo |, x <0.14 
) = (02) 


1 nt Kee Ohul 4 ee 


The relation for critical mass flux, equation (4l), is 


coupled with the two-phase momentum equation; 


CISC as PCI) 25) 28 Ge Oe Nii Uy 


O oO j 


(43) 


= c7 1x Vv + v 
O 


2 
8, - x) /2e. e 


a 

Solving equations (41) through (43) simultaneously 
generates the Henry-Fauske tables used in RELAP4. The 
tables are accessed with upstream volume stagnation pressure 
and junction stagnation enthalpy and the critical mass flux 


is returned. 
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Assumptions implicit in the use of this model are: 

(1) System flow is isentropic and one-dimensional, 

(2) The liquid phase is incompressible, 

(3) The transfer of heat and mass between phases is 
negligible, 

(4) Slip between phases is negligible and, 

(5) At a junction the vapor phase expansion can be 
described by a polytropic process. 

the Henry~Fauske model is applicable in the subcooled, 
Saturated, and superheated regions with stagnation pressure 


limits of 1 PSIA to 2400 PSIA. 


2-5 THE HOMOGENEOUS EQUILIBRIUM MODEL (HEM) (13) 


This is essentially the same as_ the Sonic model 
described in section 2.2, page 15. It differs only in how 
it is utilized by the RELAP4 code. The HEM critical 
velocity and the Sonic model critical velocity are 
calculated using the same formula, equation (6). The HFM 


critical mass flux is then calculated from, 


G = { pa} (44) 


The difference between the two models lies in how the 


density is calculated. In the Sonic model this density is 


computed from the upstream volume density by adjusting for 
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frictional losses and kinetic energy changes to the junction 
and by isentropic expansion from the Mach number just 
upstream of the junction to Mach one in the junction. In 
HEM the density is calculated by assuming that the upstream 
volume density and the pressure and enthalpy returned from 
the steam tables are stagnation properties. The junction 
density used in equation (44) is obtained by an isentropic 
expansion from the upstream volume assumed stapnation 
conditions to sonic velocity in the junction. 

The assumptions inherent in the HEM are the same as for 
the Sonic model given on page 16 with the exceptions 
mentioned in the previous paragraph. 

The HEM critical mass flux is also tabulated for use in 
RELAP& for input values of upstream volume’ stapnation 
pressures and junction stagnation enthalpy. The Homogeneous 


Equilibrium model is applicable from 1 PSIA to 3000 PSTA. 


2.6 THE MODIFIED MOMENTUM/HOMOGENEOUS EQUILIBRIUM MODEL (14) 


The Modified Momentum/Homogeneous Equilibrium model 
(MM /HEM) is the same as the HEM model except when the 


quality is less than a user specified transition quality, xe 


(the default value for this quality is X, = 0.02). 


Ti 
As mentioned in Chapter One, the simple momentum or 


inertial flow solution is quite accurate at low junction 


pressure ratios but tends to exaggerate the flow rate as 





fae! 


eritical conditions are approached. In order to compensate 
for this exaggeration the modified momentum model utilizes a 
correction to the downstream pressure based upon the 
retardation of vapor formation hy the surface tension of the 
flashing fluid. 

This more realistic downstream, or back pressure, is 


calculated from the empirical relationship, 


2 
Py z C(l - Xyp/Xp) Poat (45) 
UP 
where, 
ap SATO OOPSIA 
In the subcooled region X is set to zero and the 


UP 


critical flow rate is estimated from Bernoulli’s equation 


with the pressure drop based on the modified back pressure, 


WseeAf2g) 0 (Po 2PM : (47) 


If, however, the modified back pressure is less than or 
equal to the downstream volume static pressure then the flow 
estimate is based on inertial considerations only. In the 
transition region (0.0¢ X<X,) the critical flow estimate is 
taken as the minimum of the flow rate predicted by equation 


(47) and the flow rate calculated from, 
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1/2 

= ° 48 
"MM / HEM Cem’ (X_/Xyp? eee 
The assumptions and limitations are the same as those 


for the HEM with the additional empirical assumption of 


equation (45). 
2-7 JUNCTION FLOW RATE SOLUTION LOGIC 


If one of the critical flow models is selected for a 
particular junction the flow rate solution process for that 
junction is accomplished as follows. 

An inertial junction flow rate estimate is obtained 
from the following linearized expression derived from the 


integrated form of the momentum equation; 


Wea 7 Me t 86m - Mu?) At /1 (49) 
where all of the parameters are evaluated at the junction. 
This estimate is then compared with that for the 
particular critical flow model selected and the minimum 
identified- The fluid acceleration and junction friction 
required to force the actual calculated flow rate to the 
above identified minimum estimate is determined and_ the 
appropriate terms in the momentum equation are modified 
before simultaneous solution with the continuity and energy 


equations. 
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It is important to note that regardless of which 
critical flow model is selected, a comparison with inertial 
flow is always effected, regardless of the conditions in the 
junction and the upstream volume. Thus, if it is known 
apriori that a particular junction will not attain critical 
conditions, then selecting a critical flow model for that 
junction serves only to increase computer processinp time. 
Additionally if, for example, the Moody critical flow model 
is selected for one of the junctions connecting core 
volumes, when the fluid in those core volumes hecomes 
superheated near the end of the blowdown phase, the problem 
will abort because the Moody model is not applicable in the 


superheated regime. 
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3. OPTIMUM AREA RATIO DETERMINATION 


3.1 INTRODUCTION 


As reported in Section 1-5, the RELAP4 code in its 
present configuration, incorrectly computes stagnation 
properties. In order to alleviate this problem, the 
upstream volume flow area must be increased in order to 
prevent supersonic velocities from occurring in the volume 
immediately upstream from the break junction. 

This chapter presents the results of a study done to 
determine the optimum area ratio (AR) (ie. break junction 
flow area/ upstream volume flow area) for this problem. The 
optimum area ratio was found to be 0.7 (ie- the flow area 
for volumes V5 and V21 must be increased from 3.67 to 5.24 
square feet). In the chapters that follow, junctions, 
volumes, and heat slabs will be referred to by the letters 
J, V, and S, respectively, followed by the number of that 
junction, volume, or heat slab as defined in Figure 1 page 
7- For example, J27 refers to the break junction on the 
vessel side of the break, V5 refers to the volume 
immediately upstream of break junction J27, and V21 refers 


to the volume immediately upstream of break junction J28. 
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3.2 GENERAL 


This study made use of the Hope Creek BWR system as 
modeled for RELAP4 by EG&G Idaho, [Ince with Fvaluation Model 
controls applied (RELAP4-EM). The input deck for this 
Evaluation Model is presented in Appendix A. The critical 
flow model required for the Evaluation Model is the Henry- 
Fauske model in the subcooled regime and the Moody model in 
the saturated regime. The other Evaluation Model controls 
are irrelevant to this area ratio study and discussion of 
these controls will be delayed until Chapter Four. 

In order to understand the system transient response to 
be described in the following section, the control action 
inputs used in the problem must be explained. 

The problem starts at time 0.001 seconds when the break 
occurs. At this time the break junctions J27 (vessel side) 
and J28 (pump side) are opened. At time 0.002 seconds, J25 
(the pre-break recirculation flow junction) closes’ to 
complete the 200% break. Also at this time the feedwater 
valve, J30, closes, the two recirculation pumps, modeled by 
volumes V6 and V9, are shut off and begin to coast down, and 
a reactor SCRAM is initiated. At one second into~ the 
transient the main steam stop valve, J3l, begins to close. 
This is a slow closing valve and is not fully closed until 
four seconds into the transient. The core spray system 


(CSS), J32, and the low pressure coolant injection system 
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(LPCIS), J33, actuate 26.5 seconds after the mixture level 
in the lower downcomer, V4, reaches 21.5 feet. As a 
reference, the bottom of the lower downcomer is 7.908 feet 
below the bottom of the active core. Flow through these 
systems does not begin, however, until the volumes to which 
they are discharging (VI for the CSS and VIO for the LPCIS) 
depressurize to 304 PSIA and 310 PSIA, the rated discharge 
pressures of the CSS and LPCIS' pumps, respectively. The 
automatic depressurization system (ADS), J26, is not 
activated until 120 seconds after the mixture level in the 
lower downcomer reaches 21.38 feet and thus is not a factor 
in this study. 

In order to determine the optimum area ratio, six 
RELAP4 runs were completed with area ratios of bel. O4/9. 
0.8, 0.7, 0.6, and 0.5. The temporal response of seven key 
parameters, listed below, to the imposed LOCA transtent was 
studied; 

1) Break flow rate, vessel side of break-J2/7, 

2) Break flow rate, pump side of break-J28, 

3) Core inlet flow-J29, 

4) Lower downcomer-V4, mixture level, 

5) Lower plenum-Vl1ll, pressure, 

6) Clad surface temperature in the hottest (center) 
heat slab of the core=-S23, and 

7) Heat transfer coefficient of the fluid adjacent to 


the hot slab of the core-S23. 
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In general, the temporal parametric responses were 
bounded by the 1-0 and the 0.7 area ratio- runs. For 
purposes of clarity only these two cases will he compared 
for all of the seven parameters considered. For break flow, 
J27, and core inlet flow, J29, all six cases will be 
compared for the time interval where the greatest (and most 
important) differences occurred. 

In all of the figures presented in this chapter, the 
curves are identified by numbers (1-0, 0-9, 0.8, etc.) which 
correspond to the area ratio used to produce the data for 
that curve.e The data used to produce the figures in this, 
and subsequent chapters, was extracted from the RELAP4 
Plot/Restart tape rather than from the printed output. A 
Fortran IV program written to extract this data from the 


tape is presented in Appendix B. 


3.3 ANALYSIS OF RESULTS 


The most important parameter with respect to blowing 
down of the reactor vessel is the rate at which fluid exits 
the break. For a 200% break there are two break flows of 
Significance; the flow from the vessel side of the break 
(J27) and the flow from the pump side of the break (J28). 

Figure 3 shows the flow rate out of the vessel side of 
the break. When the break opens, the flow rate increases 


Sharply, then decreases when the flow conditions become 
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critical, or choked. As the escaping fluid increases in 
steam quality, the flow rate decreases slightly until the 
main steam Line closes at four seconds. The flow rate then 
increases momentarily- The flow rate remains’ relatively 
constant until the lower downcomer empties, allowing steam 
to escape from the break- As can be seen from Figure 4, the 
lower downcomer empties at about eleven seconds for the AR = 
0-7 case and at about 14.5 seconds for the AR = 1.0 case. 
These times correspond to the sharp decreases in break flow 
seen in Figure 3. Following the emptying of the lower 
downcomer, the flow rate decreases steadily due to the 
gradual equilization of pressure between the containment and 
the reactor vessel. 

The inset to Figure 3 reveals that the AR = 1.9 case 
initially peaks at a higher flow rate than does the AR = Q.7 
casee As will he seen, similar results were obtained for 
initial break flow on the pump” side of the break (Figure 6, 
page 42). <A satisfactory explanation for this pecultar 
result could not be found. While the difference tin this 
initial bhreak flow tis substantial, it occurs for such a 
short time that it is insignificant. 

As can be seen, the major differences between the two 
runs lie in the value of the choked flow rate attained 
immediately following break initiation. For essentially the 
same driving pressure differential the AR * 0.7 case yields 


a break flow rate of about 23,000 LBM/SEC compared to only 
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16,000 LBM/SEC for the AR = 1.0 case. This’ results in the 
lower downcomer emptying faster for the AR = 0.7 case. 
After the lower downcomer empties, the two flow rates remain 
essentially the same, with the AR = 1.0 case being slightly 
greater due to the higher junction pressure differential 
which in turn is due to a lower break flow rate in the 
initial fourteen seconds of the blowdown. 

Figure 5 shows the lower plenum pressure which ts 
essentially the driving pressure for the break flow. This 
pressure decreases slowly during the first four seconds of 
the transient because the loss of water from the break and 
the loss of steam out of the main steam line are just 
slightly greater than can be compensated for by the 
production of steam in the core. After the main steam line 
closes at four seconds the pressure increases’ slowly until 
emptying of the lower downcomer allows steam to escape from 
the break. From here on the pressure decreases’ slowly as 
the reactor vessel and containment pressures equalize. 

The differences between the two cases are due mainly to 
the time it takes for the lower downcomer to empty. Thus 
the AR = 1.0 case yields a higher maximum pressure due to 
the increased time hetween main steam line shutoff and 
emptying of the lower downcomere Following the initiation 
of steam flow out of the break the two cases’7 show 


depressurization rates which are essentially the same. 
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The flow out of the pump side of the break (J27) is 
shown in Figure 6. It should be noted that J28 is modeled 
for a positive flow into the system. For conventence, 
Figure 6 is plotted showing positive flow out of the break. 
The flow characteristics on the pump side of the break are 
significantly different from those on the vessel side of the 
break due to the presence of the recirculation pump which 
begins to shut down at 0.002 seconds- The break _ flow 
initially increases at the same rapid rate as did the vessel 
side of the break and then decreases’ rapidly. This is 
because the flow is being restricted by the recirculation 
pump and the jet pump drive nozzles- After ahout 1.5 
seconds the break flow (J28) decreases slowly and is 
oblivious to other changes in the reactor system. 

As can be’ seen, the differences between the two cases 
shown in Figure 6 are slight- The inset in Figure 6 shows 
that the peak flow rate is about 20% higher for the AR = 1.0 
casee As pointed out earlier, this difference is 
insignificant because it occurs for such a-~ short period of 
time. 

In all, it is apparent that the area ratio change had 
very little effect on the flow rate from the pump” side of 
the break. This is due to the lower break junction velocity 
caused by the flow restrictions of the recirculation pump 


and jet pump drive nozzles. 
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Figure 7 shows the core inlet flow (J29). Immediately 
after the break, the core flow begins to decrease as the 
recirculation pump on the side opposite the break coasts 
down and, as fluid normally flowing vertically upward 
through the core, is diverted to the break through the jet 
pump on the break side. This steady decrease continues 
until the jet pump intakes are uncovered. The fet pump 
intakes are located at the tops of volumes V19 and V20, and 
uncover when the mixture level in the lower downcomer 
reaches 15.7 feet.e Reference to Figure 4, page 38, reveals 
that the jet pump intakes uncover at about eight seconds for 
the AR = 0.7 case and at about ten seconds for the AR = 1.0 
case.e- These times correspond to the sharp decreases in core 
flow seen in Figure 7. Uncovering the jet pump intakes 
causes the remaining pressure differential across the core 
(due to the recirculation pump in the intact loop which is 
still coasting down) to diminish, and the core flow then 
reverses. At the time of core flow reversal, any liquid 
remaining in the core is replaced by steam from the upper 
plenum and the core dries out. The core flow remains 
negative until the water in the lower plenum begins to flash 
to steam as the system pressure drops to below saturation. 

Flashing of the water in the lower plenum causes the 
lower plenum pressure to decrease at a Slower rate than does 
the upper plenum pressure, thus the core inlet flow reverses 


again-e The core flow then decreases as the liquid inventory 
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in the reactor system decreases. The flow again turns 
negative and remains near zero for the remainder of the 
transient. 

The differences between the two caseS cAan again he 
attributed to the differences in the initial break flow 
rates discussed earlier. The lower hreak flow rate for the 
AR = 1.0 case causes the jet pump drive nozzles to uncover 
later (ten Seconds versus eight seconds), and thus the first 
core flow reversal occurs’ later. Lower plenum flashing 
starts when the system pressure drops below about 1000 PSIA, 
the saturation pressure for about 545 degrees F (the average 
water temperature in the lower plenum). Referring to Figure 
5, page 40, the system pressure drops’7 below 1000 PSIA at 
about 11.5 seconds for the AR = 0.7 case and at about 
sixteen seconds for the AR = 1.0 case. These times 
correspond to the second core reversals seen in Figure /7. 

The temporal response of clad surface temperature is 
shown in Figure 8. The data for this graph is obtained by 
taking the surface temperature of the hottest heat’ slab 
($23), which is in the center of the core. It should he 
noted that this system is modeled with the core heat slabs 
representing an average fuel assembly. Thus the surface 
temperature of heat slab $23 represents the hottest fuel 
clad temperature of an average fuel assembly. 

A more accurate core model, such as the model by 


Hendrix (9), divides the core into two regions, one 
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representing the average fuel assembly, and the other 
representing the hot fuel assembly. While this more 
sophisticated model will more accurately predict core 
temperatures, it will also increase computer processing time 
and would have added little to this analysis, which is 
concerned with relative, rather than absolute, parametric 
responses. 

Figure 8 indicates that initially the clad surface 
temperature follows’ the fluid saturation temperature 
(approximately 565 degrees F) closely until the steam 
quality in the core becomes appreciable. At the time of 
first core flow reversal the temperature of the cladding 
surface increases sharply as the core dries out and the heat 
transfer coefficient of the fluid adjacent to the heat slab 
decreases. Figure 9 shows the heat transfer coefficient 
(HTC) of the fluid adjacent to S823. The HTC decreases 
Slowly until the first core flow reversal where it drops to 
nearly zero due to dryout of the coree The time for this 
sharp decrease in the HTC corresponds to the time of first 
core flow reversal (Figure 7) and to the sharp increase in 
clad surface temperature (Figure 8). Refering back to 
Figure 8, the clad surface temperature continues its rapid 
increase until the time of lower plenum flashing and second 
core flow reversal- There it levels out somewhat due to the 
liquid flow through the core caused by lower plenum 


Flashing. Following the flow spike caused by lower plenum 
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flashing the surface temperature of the cladding again 
increases, levels out, and remains relatively constant until 
the core spray system (CSS) becomes’ effective. (The CSS 
does not become effective in the forty seconds of transient 
time shown). 

As shown by Laird (19) and Bruch (20) if the transient 
is allowed to progress beyond forty seconds, the _ clad 
surface temperature will gradually decrease until ahout 125 
seconds, when it will begin to increase. The peak claddingp 
temperatures will not be reached until well into the reflood 
phase of the accident. RELAP4-FLOOD is not available for 
BWR s and therefore the peak cladding temperature cannot he 
predicted using RELAP4 in its present configuration. 

Again, the differences between the AR = 0.7 clad 
surface temperature case and the AR #= 1.0 case lie tin the 
initial hreak flow differences which directly influence the 
time of first core flow reversal (Figure 7). As pointed 
out, the times of first core flow reversal correspond to 
drying out of the core and to the resultant sharp increases 
in cladding surface temperature seen in Figure 8. Following 
first core flow reversal the clad surface temperature for 
both cases follow nearly the same trends with the AR = 1.0 


case remaining below the AR = 0Q.7 caSe- 
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3-4 OPTIMUM AREA RATIO SELECTION 


Figure 10 shows the hreak flow from the vessel side of 
the break for each of the area ratio cases studied. In 
order to accentuate the differences between the cases, the 
transient response is only shown for the transient time from 
one to ten seconds. 

As can be seen from Figure 10, increasing the upstream 
volume flow area (ie. decreasing the area ratio) results in 
increasing the break flow rate, with the slight exception of 
the time period from five to ten seconds for the AR = 0.5, 
0.6, and 0.7 cases. While these three cases are different 
for most of the transient shown, they never differ by more 
than four percent, and for all practical purposes’ can be 
considered to be the same. 

Similar results are obtained for core inlet flow shown 
in Figure 11 for transient times from two to eleven seconds. 
The main point of interest is the time of first core flow 
reversal which occurs earlier as the area ratio is decreased 
until the area ratio reaches 0.7. Further reductions’ in 
area ratio produce no significant changes. 

As pointed out in the introduction to this chapter, the 
parametric responses studied generally fell between the AR = 
1.0 case and AR = 0.7 case. It was also noted that results 
similar to those demonstrated by Figures 10 and 1ll were 


observed for all of the parameters studied (ie. decreasing 
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the area ratio below 0.7 resulted in no significant change 
in the response of the parameter). For these reasons, an 
area ratio of 0.7 was chosen as the optimum value for this 
problem. 

It is important to note that using the RELAP4' code 
without correcting for the stagnation properties calculation 
error (ie. using AR = 1.0) yields results which are _ not 
conservative. That is, low break flow rates and low clad 


surface temperatures are predicted. 
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4. CRITICAL FLOW MODEL SENSITIVITY STUDY 


4.1 INTRODUCTION 


A sensitivity study of the critical flow models 
discussed in Chapter Two using the Standard Model (best 
estimate) configuration of the RELAP4 code (RELAP4-SM) was 
madee A comparison was also made with the Evaluation Model 
(RELAP4=-EM). It should be recalled that the critical flow 
model for the Evaluation Model is the Henry-Fauske model in 
the subcooled regime and the Moody model in the saturated 
regime. The input data decks used in this study are 
discussed in Appendix A. 

The results of the area ratio study presented in the 
previous chapter were applied to the sensitivity study runs. 
That is, the flow areas for volumes V5 and V21l1 (the volumes 
immediately upstream from the break junctions, J27 and J28, 
respectively) were increased from 3.67 to 5.24 square feet. 
The Evaluation Model run presented in this chapter is the AR 


= 0.7 case from Chapter Three. 


4.2 ANALYSIS OF RESULTS 


Figure 12 shows the break flow rate, for the first 


twenty seconds of the transient, from the vessel side of the 
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break (J27) for the six cases’7~ studied. The transient 
response for this parameter is uninteresting from twenty to 
forty seconds and is similar to the response shown in Figure 
3, page 36. In Figure 12, and in the figures and discussion 
that follows, the different cases are identified hy the 


symbology presented in Table 2. 


TABLE 2. CRITICAL FLOW MODEL SYMBOLOGY 


CRITICAL FLOW MODEL SYMBOL 
SONIC 5 
HOMOCENEOUS EQUILIBRIUM HEM 
MODIFIED MOMENTUM/HOMOGENEOUS EQUILIBRIUM MMHEM 
HENRY=-FAUSKE HF 
MOODY M 
EVALUATION MODEL EM 


Figure 12 reveals that the EM yields the highest 
average break flow rate for the initial period of the 
blowdown, and tis followed by the Moody, the Henry-Fauske, 
the MMHEM, the HEM, and the Sonic models. 

With the execption of the Moody model, the results 


depicted in Figure 12 are as would have been expected based 
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on the discussion on break flow given in Chapter Three. 
That is, following the rapid decrease in break flow 
characteristic of all of the models at about eleven seconds, 
the Moody model exhibits a peculiar jump in break flow at 
about 11.5 seconds. 

Analysis of the Moody data indicates that the lower 
downcomer empties at 11.2 seconds, allowing steam to escape 
from the break. At 11.6 seconds, the lower downcomer 
refills to a level of about 1.2 feet, and the flow out of 
the break again hecomes mostly liquid- The lower downcomer 
empties again at i2-6 seconds, and remains empty for the 
remainder of the transient. This oscillation of lower 
downcomer water level is caused by increased steam 
generation in the core starting at about 11-4 seconds. The 
increased steam generation forces water from the lower 
plenum, through the jet pump intakes, and into the lower 
downcomer-e The increase in core steam generation rate is 
caused by a return to nucleate boiling following first core 
flow reversal and core dryout. The return to nucleate 
boiling will he explained in the discussion of the heat 
transfer coefficient which follows. 

The oscillation of lower downcomer mixture level, 
described above for the Moody model, does not occur for the 
other five cases. 

The break flow from the pump side of the hreak (J28) 


showed only slight differences (less than two percent) 
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between the six cases studied. All six caseS were very 
similar to the response shown in Figure 6, pape 42, and thus 
are not presented. 

Table 3 shows the times for uncovering the jet pump 
drive nozzles and the times for emptying of the lower 
downcomer for the six cases studied. (Recall from Chapter 
Three that the jet pump drive nozzles uncover when the 
mixture level in the lower downcomer' reaches 15.7 feet, 
which corresponds to amixture level of about 8.6 feet in 
the active core). The times presented in Tahle 3 are a 
direct result of the initial break flow rates depicted in 


Figure 12. 


TABLE 3. TIME TO UNCOVER JET PUMP DRIVE NOZZLES AND 


TO EMPTY LOWER POWNCOMER 


CRITICAL FLOW TIME(SEC) TO UNCOVER TIME(SEC) TO EMPTY 
MODEL JET PUMP DRIVE NOZZLES LOWER DOWNCOMER 
EM 7-5 10.8 
M 7a 11.2 
HF 7.7 11.4 
MMHEM 8.1 12.8 
HEM 8.8 14.0 


S 9.8 14.6 
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As discussed in Chapter Three, the time that’ the jet 
pump drive nozzles uncover should correspond to the time of 
first core flow reversal, and the time that the lower 
downcomer empties should correspond to the time of second 
core flow reversal. This is indeed the case, as can he seen 
in Figure 13, which shows the core inlet flow (J29) for the 
six cases studied from zero to twenty seconds into the 
transient. The remainder of the transient was’ similar to 
the EM response shown in Figure 7, page 44 That is, the 
core flow remained near zero with small oscillations due to 
varying core steam generation rates. 

The transient responses seen in Figure 13 are as would 
have been expected from the discussion on core inlet flow 
given in Chapter Three, apain with the exception of the 
Moody model. The relatively small flow spike at the time of 
second core flow reversal for the Moody case is due mainly 
to the return to nucleate boiling discussed earlier. This 
return to nucleate boiling results in increased amounts of 
steam being generated in the core which counteracts’7 the 
effects of lower plenum flashing on core differential 
pressure. The end result is the smaller flow spike seen at 
the time of second core flow reversal in Figure 13. 

Figure 14 shows the heat’ transfer coefficient (HTC) at 
the hottest heat slab in the core for the six cases studied. 
It should be noted that, while the HTCs depicted in Figure 


14 are for the hottest heat slab in the core, for the HEM, 
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Sonic, and Moody models the hottest heat slab turned out to 
be heat slab S24 (the slab immediately above the center heat 
slab, $23) and for the other three models S23 was” the 
hottest heat slab. This result is explained as follows. 

One of the main differences between the Fvaluation 
Model (RELAP4~-EM) and the Standard Model (RELAP4~SMN) is in 
the heat transfer logic used in the two cases. For the case 
of RELAP4{~EM, it is stated in Appendix K to 190 CFR 50 that 
"after critical heat flux is first predicted at an axial 
fuel rod location during blowdown, the calculation shall not 
use nucleate boiling heat transfer correlations at that 
location subsequently during the blowdown even if the 
calculated local fluid and surface conditions would 
apparently justify the re-establishment of nucleate 
boiling". RELAP4&-SM, = on the other hand, allows re-~ 
establishment of nucleate boiling if the conditions warrant. 
For the HEM, Sonic, and Moody model cases, at the time of 
first core flow reversal, the conditions were such that 
nucleate boiling was re~established at heat slab $23, while 
film boiling continued at heat slab S24. Nucleate boiling 
conditions persisted for about ten seconds at heat slab $23 
for both the HEM and Sonic models and for about’ three 
seconds for the Moody model. This resulted in the 
temperature of heat slab S23 being lower than that for heat 
slab S824 in all three cases for the remainder of the 


transient. 
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Figure 15 shows the cladding surface temperature at the 
hottest heat slab in the core for the six cases studied. 
Again, for the HEM, Sonic, and Moody model cases, the 
hottest heat slab was S24, while for the other three cases 
S23 was the hottest. 

Table 4 shows the peak cladding surface temperature 


attained during the blowdown for each of the six cases 


studied. 
TABLE 4. PEAK CLADDING SURFACE TEMPERATURE 
CRITICAL FLOW PEAK CLADDING SURFACE 
MODEL TEMPERATURE (DEG F) 
EM 1349 
HF 919 
M 909 
MMHEM 898 
HEM 850 
S 846 


As can be’ seen from Figure 15 the surface temperature 
decreases rapidly for the Moody model case at ahout eleven 


seconds. This decrease corresponds’ to the re-establishment 
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of nucleate boiling which, for this model only, occurs in 
heat slab S24, as well as in heat slab S23. This return to 
nucleate boiling at heat slab S24 is characterized hy the 
slight increase in tthe HTC seen at about eleven seconds in 
Figure 14. With the exception of this one peculiarity for 
the Moody model, the surface temperature transient responses 
for all six cases follow trends that are consistent with the 


discussion of clad surface temperature in Chapter Three. 
4.3 CONCLUSIONS 


From the analysis of the previous section it is 
apparent that the Evaluation Model, as expected, produced 
the most conservative results. As pointed out, this is 
primarily due to the different heat transfer logic in the 
RELAP4-EM and RELAP4=-SM configurations. 

Based on initial break flow rate (Figure 12, page 55) 
it is clear that the Moody model is most conservative, 
Peiiiceea by Henry-Fauske, MMHEM, HEM, and Sonic. Note, 
however, that the Moody and Henry-Fauske models differ by 
only about five percent. 

Based on peak cladding surface temperature attained 
(Table 4, page 63) one would conclude that the Henry-Fauske 
model is more conservative than the Moody model, followed by 
the other models in the same relative order as ahove. Note 
again, that the difference between the Moody and Henry- 


Fauske models is slight (only about one percent). 
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For these reasons, the Moody and Henry~Fauske models 
are rated as equivalent hased on conservative results, 
followed, in order of their relative degrees of 
conservatism, by the MMHEM, the HEM, and the Sonic critical 


flow model. 
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3. SUMMARY 


Dei RESULTS 


This study utilizes the RELAP4/MOD5 Computer Code to 
analyze the effects of critical flow modeling on the 
thermal-hydraulic transient response of a General Flectric 
boiling water reactor to a major primary coolant line 
rupturee- Included in the study is a presentation of the 
equations, assumptions, and limitations of the critical flow 
models available for use in RELAP4. Additionally, an 
evaluation of a temporary solution to a RELAP4 coding error 
associated with stagnation properties calculation is 
presented. 

The results of this investigation indicate that: 

(1) A solution to the stagnation properties 
calculational error that exists in the RELAP4 code, when 
applied to the Evaluation Model, (RELAP4-EM) provides a 
conservative evaluation, relative to the Standard Model 
(RELAP4-SM), regardless of which of the five critical flow 
models available in the code is selected; 

(2) Of the five critical flow models available for use 
in RELAP4, the Moody model and the Henry-Fauske model are 
nearly equivalent in their relative deprees of conservatism. 
Following these models, in order of relative degree of 


conservatism, are the Modified Momentum/Homogeneous 
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Equilibrium model, the Homogeneous’ Equilibrium model, and 
the Sonic model; 

(3) Ta order to best alleviate the stapnation 
properties coding error, the flow areas for the volumes 
immediately upstream from the recirculation line rupture 
should be increased such that the area ratio (break flow 
area / upstream volume flow area) is 0.7. 

All of the critical flow models in RELAP4, with the 
exception of the Sonic model, require the use of stagnation 
properties for estimation of flow rate under critical or 
choked conditions. The RELAP4 code, in its present form, 
incorrectly computes stagnation properties when the flow 
rate is such that the code predicts supersonic velocities in 
the volume immediately upstream of the junction where 
critical flow conditions exist. The results of this 
investigation indicate that the optimum solutton to this 
problem is to arbitrarily increase the upstream volume flow 
area such that the area ratio is 0.7. 

The Hope Creek boiling water reactor, included as a 
sample problem in the RELAP4 Users Manual, was used _ to 
analyze the effects of critical flow modeling on the 
predicted blowdown response of this BWR to a double-ended 
guillotine rupture of the primary coolant’ recirculation 
line- Applying the above mentioned stagnation properties 
calculation problem "fix" to RELAP4-EM provides a 


conservative evaluation relative to RELAP4=-SM, regardless of 
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which of the five critical flow models available in the code 
is selected. Applying the "fix" to the Standard Model tt 
was found that the Moody and Henry-Fauske critical flow 
models provide nearly equivalent evaluations in terms of 
peak cladding temperatures attained during the bhlowdown. 
The other three critical flow models, in order of their 
relative degrees of conservatism, were found to be the 
Modified Momentum/Homogeneous Equilibrium model, the 


Homogeneous Equilibrium model, and the Sonic model. 


5-2 RECOMMENDATIONS FOR FURTHER STUDY 


Based on this investigation the following recommenda- 
tions for future research are presented: 

(1) While this investigation rated the critical flow 
models by their relative degree of conservatism, comparison 
of the results of RELAP4 analysis with LOFT, SFMISCALE, or 
other large scale experimental data would provide valuable 
insight as to which of the models most accurately predict 
flow rate under critical or choked conditions. 

(2) The heat transfer logic employed by the RELAP4 code 
appears to have a very significant effect on the overall 
performance of the code. Analysis of the effects’ of 
changing heat transfer and critical heat flux correlations 


would be invaluable to future users of the code. 
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APPENDIX A 


INPUT DATA DECKS 


This appendix provides a listing of the input data 
decks for the Evaluation Model (RELAP4-EM) and the Standard 
Model (RELAP4-SM) used in this investigation. The data 
consists of title cards, data cards, and comment cards. The 
comment cards are prefixed with an asterisk (*). 

Table 5 is the input deck for RELAP4-EM that was used 
in the area ratio study presented in Chapter Three. The 
only two cards that were changed to produce the different 
area ratio runs were the volume data cards 050051 = and 
050211. Card 050051 is is the data card for volume V5, and 
ecard 050211 is the data card for volume V21. The underlined 
number on these cards is the volume flow area. This 
quantity was varied from 3-67 square feet for the area ratto 
(AR) = 1.0 case to 7.34 square feet for the AR = 0.5 case. 
The example shown in Table 5 is the AR = 0.7 case. 

Table 6 is the input deck for RELAP4-SM used in the 
sensitivity study of Chapter Four- Note that the volume 
flow areas for V5 and V21 (cards 050051 and 050211) are 5.24 
square feet, the optimum value found in Chapter Three. The 
only three cards changed in this deck to produce’ the 
different critical flow model runs were the title card (the 
first card in Table 6), and the junction data cards 080271 
and 080281, for break junctions J27 and J28, respectively. 


The underlined quantities on these cards are _ the critical 
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flow model selection indices ICHOKE and JCHOKE. JCHOKE is 
the first underlined quantity, ICHOKE is the second. The 
values of JCHOKE and ICKOKE determine which critical flow 
model is to be used for that junctione Table V, page I-14/7 
of the RELAP4 Users Manual is used to determine the values 
of JCHOKE and ICHOKE for the desired critical flow model. 
In the example shown, for junctions J27 and J28, the values 
are JCHOKE = 2, ICHOKE = | which indicate that the Sonic 
critical flow model is to be used for junctions J27 and J2®&. 
Note that all other junctions are set up to use the inertial 


flow model (ie. JCHOKE = -1, ICHOKE = 1). 
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TABLE 5. EVALUATION MODEL INPUT DATA DECK 


=BWR EVALUATION MODEL 200% RBREAK USINC 7 X 7 FUEL 
*234567890123456789012 3456 789012 345678901234567890 

*10001 LNDOM EDI NTS TRP VOL BUB TDV JUN PMP CFV NLK FLL SL&AV GOM MAT COR 
010001 -2 ) 5 10° 23 360 33°2°46 2 4 25310°5 5-0 1 

010002 3388. 1.9 

010003 0 l 1 l 

020000 AP 11 JW 27 JW 28 JW 29 JW 16 ML 4&4 SR 23 CR 23 TD 27 


050181 1045.86 532.67 -l- 534.53 12.35 12.35 0 63.635 0. 17.92 1 


030010 1 100 ®, 0 -001 -00001 -l 2000. 
030020 5 8 3 0 01 000001 l. 
030030 20 10 2 0 01 -000001 60. 
030040 20 25 l 0 ~01 -000001 130. 
030050 20 25 l 0 Ol -00001 2000. 
030003 50 
040010 1 1 0 0 60.0 O. 
040020 1 -4 3 0 40.0 0. 
040030 2 1 0 0 0.0 0. 
040040 3 1 0 0 0.001 QO. 
940050 4 1 0 0 002 0. 
040060 5 1 0 0 1.0 0. 
040070 6 1 0 0 1.E6 0. 
040080 7 =5 4 0 213182065 
040090 8 1 0 0 -001 1.E6 
040100 9 =-5 4 0 21.38 120. 
RSRXXX VOLUME DATA CARDS 
050011 2 0 1039. =I. <O91632 1548. 22.625 22.625 0 226. 0. (29.92 0 
050021 2 0 1028.28 -1. .002448 6703. 12.25 5.27 0 497. 0. 41.75 0 
050031 1 0 1024.98 -l. .999 3777. 21.25 21.25 0 497. 0. 51.25 O 
050041 3 0 1031.93 532.3 -l. SU7?%« $2023 32023 0.165. 0. 10.02 0 
050051 0 0 1024. 532.3 -l. 125025 208 26. 9 5.24 0. “14. 9O 
050061 0 O 1126-17 533.25 -l. 63. 36 3. 04.90 0. =-14. O 
050071 0 0 1226.17 533.25 -le- 366. 40.53 460.53 0 3.2 OO. =14. O 
050081 0 0 1024. 5$2.3)9 w-l. 136. 26. 26%. 0 3.67 0. -14. 9 
050091 0 0 1126.17 533.25 -l. 63. 3. 3. 04.0 0. =-14 £49 
050101 90 O 1226.17 533.25 =-l. 366. 40.53 40.53 0 3.2 O. =-14 9 
050111 2 O 1062.16 532.53 -=l-e 2131-5 17-20 17.20 O 120. OQ. 0. 23 
050121 2 O 1054.67 549.337 -1. 223. 2.85 2.85 0 81.091 .0473 17.82 13 
050131 2 0 1052.97 -1. 2.037157 223. 2-75 2.75 0 81.091 .0473 20.67 14 
050141 2 90 1051.52 -1. .050567 81.09 1. l. 0 81.091 .0473 23.42 15 
050151 2 0 1049.87 -1. .091994 223. 2.75 2.75 0 81.091 .0473 24.42 16 
050161 2 0 1047.65 wl. 139875 223. 2-75 2.75 0 81.091 .0473 27.17 1 
050171 2 0 1054.28 532.53 -l. 942.47 18-12 18.12 9 52.301 0. O. 18 
2 0 
050191 0 0 1059.2 532.53 <1. 115. 15.7 15.7 O 19.69 OO. 10.02 0 

050201 0 0 1059.2 532.53 -1l. 115. 15.7 15.7 0 19.69 O. 10.02 0O 

050211 0 0 1030.8 532.30 -l.- 10-75 2.22 Paz 0.5.24 O. 10. 0 

050221 4 0 14.7 100. 0.6 3.43E6 250. 0.0 0 1000. 100. -10. 0O 

050231 2 0 1055.48 532.53 -l. 20.5 Le25 1625 0 861.09 OO. 16.77 12 


*6000 LIQUID LEVEL CALCULATION CARD 
VeE0000 Tb 23 12°13 14 15 16 17 18 19 20 1 
060001 0.0 0.0 0.0 0.0 0.0 
*60002 WALLIS SLIP CORRELATION CARD 
*60002 WALSC1l WALSC2 
*6XXX1 BUBBLE DATA CARDS 
*6XXXL ALPH VBUB 
060011 1.0 0.0 
060021 0.0 3.0 
060031 0.8 3.0 
060041 0. 1.£E6 
*60001 SLIP VELOCITY CARD 
*BXXXY JUNCTION DATA CARD 
*2346567890123465678901234567890123456789012 3456789012 34567890123456789012 


080011 1 #2 0 0 29583. 191.1 52-17 0. Of. O- 10300.0 0. 11:0 1-0 0 
080021 2 3 00 4074.3135 312. 52.25 0. 0. O- 1030 0.0 0. 1100.0 0 
080031 2 40 0 25444. 165. 42.0 0. 0. 0- 1903 00-0 0. 1100.0 0 
080041 4 20 0 0 10041.5 3.112 25. 30. .2084 1-17 2 002 0.0 0. 1100. 0 
080051 4 21 0 0 4750. 3.67 11.0 0. O. O- 10300.0 0. 11 00.0 0 
080061 5 6 =1 0 46750. 3.67 -13. 0. O. 0O- 00 3 00.0 0. 1190 0. 0 
080071 6 7 1 0 4750. 3.2 -13. 0. 0. Q. 0 0 3 0 0.0 Q. ll 0. 0. 0 
080081 7 20 00 4750. -538 25. 47. 2.2373 6.8 0002 0. 0. 110 0. 0 
080091 4 8 00 4750. 3.67 11.0 0. O-. O- 103 00.0 0. 1100. QO 
080101 8 9% =2 0 4750. 3.67 -13. O-. O-. 0O- 00300. QO. 11:00. 9 
080111 9 102 0 4750. 3.2 <-13. 0s 06. 0. 0-0 3 OD. OO. Il 0 0. 6 
080121 10 19 0 O 4750. .538 25. 47. .2373 6.8 00020. 0. 1100. 0 
080131 4 19 0 0 10041.5 3.112 25. 30. .2084 1.17 2 002 0. Of. 11:0. «0. 9 
080141 23 12 0 0 26639. 81.09 17.92 0 0. 0.003090. 0. 11:0 1.00 
080151 12 13 0 O 26639. 81.09 20.67 0. 0. 0.00300. 0. 110 1.0 0 
080161 13 14 0 0 26639. 81.09 23.42 0. 0. 0. 090300. 0. 11:0 1.0 0 





Bs 


TABLE 5- (cont. ) 
080171 14 15 0 0 26639. Cle09 242.42 06. 0. 0. 0 G0 3 0 0. 0. Il 0 1.0 0 
VBOLTEL [5 16 0 0 26639. 81.09 27.17 0. 0. 0. 0 03 0 0. O« 11 0 1.0 0 
O80191 16 190 0 26639. 37.13 29.92 0. 0. 0.00300. 0. 11:0 1.0 0 
080201 11 17 0 0 2944. 2-246 17.0 QO. 0 02 190<30t1.0 0. 11 0 0-0 0 
080211 17 18 0 0 2944. 26-20 18.02 0. 0. 0. 00300. 0. i1 0 0.0 0 
080221 18 1900 2944. 43.635 30.02 0. 0. 0. 00300. 0. Ll O 1.0 0 
080231 19 11 0 0 14791.5 19.69 10.12 0. O. O 1909300. 0. 11:0 0.0 0 
080241 20 11 00 14791.5 19.69 10.12 0. O«. 0. 10300. 06. 11:0 90.0 0 
OBOZ5l 21 5 O 1 4750. Bienes 11.0 O. O. 0 09030 0. O- 11:0 0.0 9 
O82026l 3 22 O33 G. «484 55-0 IL. le le 000090. O- LL O9N.0 O 
080271 21 22 04 O. 30 67 Il. OO. Ils «5 0 0 200. .9 Ii 0 0.0 6 
080281 22 5 04 0. 3.67 Il. Osc «<3 Ie 0 0 2 0 0. «9 11 0 0-0 0 
O80Z29t L£1 23 0°70 26639. 20.32 L720 0.00705. ) 0 3 0 0. 0. 11 0 1.0 0 
080301 0 4 1 0 4074.3135 4.712 40. 0. 0. 0. 00300. 0. Ill 0 0.0 0 
080311 0 3 2 2 =4074.3135 14.75 2- O. O02 O<f. 0 03 0 0. 0. 11 00.0 0 
Oe0321 0 13.0.0. 1.0 30. OO. O. 0. 00300. 0. 1190 0.0 0 
080331 0 19 4 0 0. 1.00 256 ls Oe 00.00-0 0:0 0.70. 11.0 0.0 0 
082000 1 
*9OXXY PUMP DESCRIPTION DATA CARDS 
090011 3 4 OQ 1 0 1668. 1.0 44943. 710.3 22200. 20000. 47.17 Oo. OQ 
090021 3 4 0 1 0 1668 1.0 44943. 710.3 22200. 20000. 47.17 0. 9O 
O91L00!1 -1l QO. OQ. ol 0. °15 .05 .24 «2.8 23 °-96 44 -98 .6 ao 
091002 -8 .9 9 8 96 «.5 1.0 OQ. 
092001 -2 0. 0. L. OQ. 
Cgs0tl 0.  O. 0. 
095021 0. QO. QQ. 
100000 0 0 16 90 
103011 if 53 0.0 1.31 Oazs Paz3 Os 1.20 Oat5 1.12 1.0 1.0 
fOsgzt tf 2 5 0.0 -0.9 225 -.5 25 -.05 cia °4 1.0 1.0 
Os0st to 73-3) =1.0 221 -.75 1.775 -0.5 1.54 -.25 1.4 Q. 1.31 
103041 1 4 5 #-1.0 2.1 -.75 1.50 -0.5 1.13 -.25 92 0. ~825 
Heso5i 1 #5 5 0. 74 e295 © 53 5 ©625 ~75 - 80 1.0 1.90 
Pe3061 1° 6 5 O. 0825 225 -79 25 815 275 -90 1.0 1.0 
103071 lL 7 5 #=1-90 -1.8 ~.75 -l. -.5 -.4 -.25 el 0. 4 
103081 l 8 5 -1.0 -1.8 -.75 -1.72 -.5 -1.54 -.25 -1.27 0. -.9 
PO309LF 2 %I 5 O. 255 025 ~63 25 71 075 - 86 1.0 1.90 
103101 2 2 5 0. -.55 o2 -2175 4 «0 Pe | - 58 1.0 1.0 
POSIT 2 3 S =1.0 tsi? -.75 be23 -.5 83 -.25 62 0. °55 
103121 2 4 5 #1.0 1.77 -.75 1.52 -.5 1.32 -.25 Daal 3 0. 1.0 
POst3i 2 5 5 0. -.8 0.2 -.375 0.5 -.l ©75 «15 1.0 »43 
POS1Gl 2° 6 °5 0. 1.0 225 - 88 25 275 ay Be 61 1.0 ~ 43 
bOsTSi>> 2.7 °5 =I. -3.45 <-.75 -2.75 -.50 <=-2.0 -.25 -1.375 0O. -.8 
103161 2 8 5 el. -3.45 <-.75 -2.62 --50 e-1.583 -.25 “1.15 0. -.955 
*ILLXXXO VALVE DATA CARDS 
110010 4 QO. 900. QO. QO. OQ. 
PEOOZO> “5 2 0 -025. 0. O«. °04 
110030 -9 0. 0 O. QO. OQ. OQ. 
110040 -+3 0.9 0. 0. 0. 0. 
*I2XXYY LEAK TABLE DATA CARDS 
120L01l1 -+2 9 30. Ol. 1000. it. 
20201 6 5 1467 O«. 1s 2425 Le 265 «667 2575 «333 3. 0. 10000. 0. 
*I13XXYY FILL TABLE DNATA DARDS 
eoxOO TRIP TYPE POINTS ICALC UNITS PORT HORX 
130100 2 1 4 1 LBS/SEC 1050. 401.11324 
130101 O. 864.6676 1. 864.6676 4.0 0. 10000. OQ. 
130200 2 2 4 5) LBS/SEC. 548. 0999 
130201 O. Of. 940. O. 1025. -276-22246 3000. -276.2246 
130300 aL 2 Ll 3 GAL/MIN 212. 0.0 
130301 0. 15150. 59. 14200. 137. 12500. 178 11550. 237. 9470. 
130302 268. 7576. 280. 5682. 287. 3788. 297. 1894. 304. 0. 1LNOQON. O. 
130400 Ta 2 9 3. GAL/MIN 202 0. 
130401 17.6 4355. 23.7 4309. 26.1 4291. 41.2 4175. 42.4 4166. 
130402 98. 3708. 102. 3672. 310. 0. 3000. OQ. 
*14000 KINFTICS CONSTANTS 
*14000 PWR MUL BIL REC UDUF 
140000 3 116. 0. l. 
*141001 SCRAM TABLE DATA CARDS 
141001 -8 4 0. 0. 0.9 O- 1.075 -.4 1.6 -4. 2-7 -10. 5-7 -28.1 8. -30. 10. =-30. 
*142001 DENSITY REACTIVITY DATA CARDS 
142001 -10 .2 11.55 -3 -8.830 .4 -6.690 .5 <-5.0l °6 -3-64 .7 =-2.51 
14209002 8 -1.55 09 -0.73 l. 0. 1.1 O. 
*143001 DOPPLER TABLE DATA CARDS 
143N01 -6 QO. 1.85 1000. 0. 2000. -1.71 3000. -3. 4000. -4.28 S000. -5.43 





TABLE 5. (cont.) 


*140XX0 REACTIVITY COEFFICIENT DATA CARDS 

140010 .148289 .148289 0. O. 

140020 .341782 .341782 OO. O. 

140030 .019858 .O19R58 O. O. 

140040 .341782 .341782 O. QO. 

140050 .148289 .148289 0. O. 

*L5XXX HEAT SLAB DATA 

DoOOrh ay 3) 20 0 Tl 11 GO. G87. 4622. 0. 0. 0. 0. 0. 0. 
150012 0. 9. 

150021 0 2 290011 LI 0. 1368. 698. O. 0. Of. O- O- OQ. 
150022 0. O. 

150031 09 42 00 11 11 0. 1916. 9772. Oe Oe Of Oe Oe Oe 
150032 0. QO. 

150041 9 11 20011 11 0. 687. 385. OO. O. O. O- OO. QO. 
150042 0. 0. 

150051 4 18 93 0 O LL Ll 1355. 1355. 239. O. Of Of O- O- OD. 
150052 0. O. 

150061 1 240 0 11 LL 2288. 2288. 286. 0. 0. O- OO. O- 0. 
150062 0. 10.42 

150071 11 23 50011 II 197- 197 65. Of Of. Oj. Ow} O- OW 
150072 0. 0. 

150081 16 1 5 00 11 IL 94. 94. gl. Of. 0. DU. 0. O. 0. 
150082 0. 0. 

150091 Il 17 6 00 11 Ll 6884. 6884. 95. 0. O-. O./. O- DO O. 
150092 OO. IL7.l 

I5010L 0 11 700 iLL IL O. 736. 46. 0. Of Oe Oe Oe Oo 
150102 0. O. 

150111 2 3 8 0 O Ll IL 4380. 4380. 184. O. Of Of Of O- O. 
150112 0. QO. 


150821 05 9 0 0 LI 11 0 450. 36. OO. Of O- OO. ODO. OD. 
50122. 0. 0. 
150131 08 90011 LI O. 450. 36-4 QO. O01. O- O- O- OO. 
150132 0. 0. 
150141 0 790011 11 O.- 838. 67. QO. Of O- Of ODO. OD] 


150142 0. 0. 

1508s 010 VCO 0 Tl Ll 0. 838. 67. O. O. O«. O. DO. A. 

150152 0. 9. 

P501Gl 12 18 10 00 Il 11 1096.8 1096.8 25.67 0. 0. 0. 0. 16375 1.375 
150162 0. 2.74 

150171 13 18 10 L O Il IL 1096.8 1096.8 25.67 0. Ow. O«/. DO] 4.125 4.125 
150172 2.76 5.49 

150181 14 18 10 I O LL Il 398.84 398.84 9.334 0. 0. O. O- 6.0 6-0 

150182 5.51 6.49 

150191 15 18 LO L OLL IL 1096.8 1096.8 25.67 0. O« Of O. 76875 7.875 
150192 6.51 9.24 

150201 16 18 10 1 O LL LL 1096.8 1096.8 25.67 0. O. O- OO. 10.625 10.625 
150202 9.26 11.99 

f>02ti 0 12 10 0 Ll Ll ©. 15175.52 178.039 0. .0468 0. «0568 1.375 1.375 
150212 0. O. 

150221 90 13 Lt tL O Ll IL Of 15175-52 178-039 O«/ .048 0. .0588 4.125 4.125 
[50222 0. 0. 

150231 0 14 1 2 0 LL LL Of 5518.372 64.741 0 .048 0. .0588 6.0 6.0 
15023920. 0. 

150241 0 15 L LO 11 IL Of 15175.52 178.039 0. .048 0. .0588 7.875 7.875 
150242 0. 0. 

150251 0 16 Lt 10411 IL O. 15175.52 178.039 0. .048 0. .0588 10.625 10.625 
150252 0. 0. 

*23456789012345678901234567890123456789012345678901234567890 

*16000 CORE SECTION DATA 


*16000 SLB NODE CLAD QFRAC 
L60010 21 79 15 0. eL7721 
[60020 22 79 15 QO. -26902 
L6n0g0 23 4 9 I5 0. 10754 
160040 24 79 15 O- 26902 
160050 25 7°9 15 OQ. 217721 
160015 0 O 96.07 .8 06 © 965 -95 00 .05 0 0 0 .08333 .9714 .01988 
160025 0 900. 8 °6 °965 «95.0 0 2.05 © 0,0 206337 0. -01988 
160035 0 900. 8 -6 965 -95 00 .05 0 00 .08333 0. -01988 
160045 0 00. 8 6 965 °95 00 .05 00 0 .08333 0. -01988 
160055 0 Q .1259 .8 26 © 965 °-95 0 0 .05 00 0 .08333 .00127 .01988 


ooo0ce 
oooc]e 





Jan 


TABLE 5. (cont.) 
*17XXYY SLAB GEOMETRY DATA CARDS 
[7Ol0t 2 3 1 6 0.0 -020071 1.0 
170102 1 23.2 -000721 0. 
170103 O 3 8 -00267 0. 
1jozol 1 2-4 3 0.0 +4 OQ: 
170202 0 5 1 -0208 0. 
170360 1} 1 5S 1 0. ee On 
170401 41 #5 #1 4. 2125 0. 
170501 1 1 5 10. a3 On 
170601 1 1 5 10. -0138 O. 
P7O7Oletl “f “5 1 0. 0625 Oc 
l7osel 1 ft 5S 1 0. 042 0. 
170901 1 1 5 10. -08 Ge 
171001 24 1 5 10. -0234 0. 
180100 -20 500. 3.341 650. 2.971 800. 2.677 950% P2499 LI00G) 72242 
180101 1250. 2.078 1400. 1-940 1550. 1-823 1700. 1.724 1850. 1.639 
180102 7000. » 1.568 2150. 1.507°° 2300. 1.457 2450. 1.415 2600. 1.382 
180103 3100. 91.323 3600. “1.333 61700: 16406 4600. 1.538 Sioo. 1.730 
180200 2 37. 41562 5400. 41562 
180300 18 32; 7-812 212% 7.992 392. 8.208 
180301 S72. 8.784 S96 9.540 032% 10.404 
180302 Ty ea 11.268 1292. 12.492 1472. 135176 
180303 1652. 13.968 1 832% 14.796 2012. 16.128 
180304 2192. 17.784 js A des 19.656 2552. 21.780 
180305 TIS26 24.048 3092. 28.908 3360. 372128 
P8@401-=5) 32-30. 212. 29.5 392. 28.3 S72. 2656 752. 24.7 
180501 -2 200. 8.33 1200. 12.92 
190100 16 32. 34.45 122. 38.35 212. 40.95 
190101 392. 43.55 752. 468 2012. 51-35 
190102 2732. 52.65 3092. 56.55 3452. 63.05 
190103 1c rage bare.) 4352. 89.7 4532. 94.25 
190104 4712. 98.15 4892. 100.1 51442 10164 
190105 8000. 1LOL.4 
190200 2 van -000075 5400. -000075 
190300 5 0.0 28.392 
190301 1480.3 34.476 1675.0 85.176 
190302 1787.5 34.476 3500.0 34.476 
190401 =7 130. 56-9 350. 60-8 450. 6263 530. 6562 620. 67-2 710. 70-2 800. 77.5 
190501 -10 68. 52-8 200. 56.7 400. 61-6 600. 64. 800. 66. 1LOOO. .67 
190502 1200. 68-4 1400. 71-8 1600. 75-8 1800. 80.6 
200101 -2 0. 3.718E-6 5000. 1.2653E-5 
200201 -2 0. 0. 5000. 0. 
200301 -4 OQ. 3-094E-6 1652. 4-.706E-6 1653. 5.389E-6 5000. 5.389E-6 
200401 -2 O- O. 5000. 0. 
200501 -2 O. O. 5000. 0. 
250001 1 2 3 0 .0615 
260001 1234666611 1LE7 11000000000 
270001 -14 2500. 0. 2200. 98.7 1620. 197. 1730. 395. 1660. 592. 1600. 790. 
270002 1540. 987. 1480. 1184. 1440. 1382. 1400. 1579. 1370. 1777. 1335. 1974. 
270003 1310. 2172. 1280. 2369. 
280001 -4 34. O. 34. 4150. 21. 200. 21. 1000. 
290001 -14 2500. 0. 2200. 98.7 1620. 197. 1730. 395. 1660. 592. 1600. 790. 
290002 1540. 987. 1480. 1184. 1440. 1382. 1400. 1579. 1370. 1777. 1335. 1974. 
290003 1310. 2172. 1280. 2369.- 
300001 -14 131. O- 131. 98-7 626 197. 47-6 395. 69- 592. 95. 790. 116. 987. 
300002 Voice HI A46- 037s 013826 19 1509s 12 ls 1777 10G. 1974.5 766. 2172. 
300003 58. 2369. 





TABLE 6. STANDARD MODEL INPUT DATA DECK 


=™BWUR BEST ESTIMATE 2002RRK CRITICAL FLOW = SONIC 
*2345678901234567890123456789012 345678901234567890 
*10001 LDM EDI NTS TRP VOL BUB TDV JUN PNP CKV NLK FLL SLBV GOM MAT COR 


010001 -2 9 5 or 23 14 °0 3924 2-4 25 (10 5 5 0.0 

010002 3388. 1.0 

020000 AP ll JW 27 JW 28 JW 29 JW 16 ML & SR 23 CR 23 TD 27 

030010 1 100 5 0 -001 -00001 sl. 2000. 

030020 5 8 70 s01 -000001 ls 

030030 20 10 2 0 Ol -000001 60. 

030040 20 25 1 0 01 -000001 130. 

030050 20 25 1 0 01 -00001 2000. 

030003 50 

040010 1 1 0 Oo 60.0 0. 

040020 1 -4 3 0 40.0 oye 

040030 2 1 0 0 0.0 0. 

0400460 3 1 0 0 0.001 0: 

040050 4 1 0 0 002 0. 

040060 5 1 0 0 1:50 0. 

040070 6 1 0 oO 1.E6 Os 

040080 7 <5 4 0 21-5 26-5 

040090 8 1 0 Oo 001 1-E6 

040100 9 <5 4 0 21.38 120. 

*5XXXX VOLUME DATA CARDS 

O50011 2 0 1039. eto 6091632 15468. 22.25°22.25 0 224%. 0... 20.92 0 

050021 2 09 1028.28 -1. .002448 6703. 12-25 5.627 0 497. O.« 41.75 0 

050031 1 0 1024.98 -1. .999 M777 21.275 21225 0 497.°0. ~$1425.0 

050041 3 O 1031.93 532.3 <1. 3177. 32-23 32-23 0 165. O. 10.02 O 

050051 0 0 1024. §32.3 <1. 125525532 6. 26. O°5.24 0. =—1450 0 

050061 0 0 1126.17 533.25 <1. 63. 3 cm 0420) Osu =1462 0 0 

050071 0 O 1226-17 533-25 -l. 366. 40.53 40-53 03.2 Of -14 0 

050081 0 0 1024. §32.30 -1. 136. 26. 26. 03667 0. =<bés “0 

050091 8) 0 1126.17 8362.5 -l. 63. 3e 3 6 6) 4.0 0. -14. 0 

OSO10L O O 1226.17 533.25 -l. 366. 40.53 40553 0 3.2 0. “=14.. 0 

OSOLIL 2 O 1062.16 532.53 -1. 2131.5 17-20 17-20 0 120. O. 0. 23 

OSOl21 2 01054267 549.3397 =I 223. 2085 2.85 0 81.091 .0473 17.82 13 

050131 2 0 1052.97 -1. .037157 223. 2-75 2.75 O 81.091 .0473 20.67 14 

050141 2 0 1051.52 -1. .050567 81.09 lL. it 0 81.091 .0473 23-42 15 

050151 2 O 1049.87 <1. .091994 223. 2.75 2-75 O 81.091 .0473 24.42 16 

050161 2 O 1047.65 -1. .139875 223. 2-75 2-75 U0 81.091 .0473 27.17 1 

050171 2 0 1954.28 532.53 <1. 942.47 18.12 18.12 0 52.301 0. O. 18 

O50181 2 0 1045.86 532-67 -1le 534653 12-35 12-35 0 43.635 O- 17.92 1 

050191 0 0 1059.2 532-53 <1. 115. 15-7 15-7 0O 19.69 0. 10.02 0 

050201 9 O 1059.2 532.53 <1. 115. 15.7 1567 0°19.69 0. 10.020 

0502010 0° 1030.8 532.30 <1. 10.75 “2.2 252 0 5.24 oO. 10. 0 

050221 4 0 14.7 100- 0.6 3-43E6 250. 0.0 1 1000. 100. -10. 0 

050231 2 0 1055-48 532.53 <1. 20.5 1625. 2225 te 81.090005 16277. 12 
*6000 LIQUID LEVEL CALCULATION CARD 

060000 11 23 12 13 14 15 16 17 18 19 20 1 

060001 0.0 0.0 0.0 0.0 0.0 

*60002 WALLIS SLIP CORRELATION CARD 

*60002 WALSCL WALSC2 

*6XXX1 RUBBLE DATA CARDS 

*6XXX1 ALPH VBUB 

060011 1.0 0.0 

060021 0.0 3.0 

060031 0.8 3.0 

060041 0. 1.E6 

*60001 SLIP VELOCITY CARD | 

*8YXXY JUNCTION DATA CARD 


*234567890123456789012 34567890123456789012345678901234567890123456789012 


So 


080011 1 #200 29583. 191-1 52-17 9. Of. O- 1-1 3 O0O-.0 O- 1O01.0 0 
080021 2 30 0 4074.3135 312. 52.25 0. O- O- Ll -1 3 00-0 0. 1 0 0.0 0 
080031 2 40 0 25444. 165-4. 42.0 QO. O- O 1-1 3 O 0.0 O- 10 0.0 0 
080041 4 20 0 0 10041.5 3-112 25. 30. .2084 1.17 2 =-1 0 2 0.0 0. 1 0 0. 
080051 4 2100 4750. 3-67 11.0 O. O. O- 1 -1 3 00.0 0O- 1900-0 0 
080061 5 6 -1 0 4750. 3.67 -13. 0. Oj. 0. O -1 3 0 0.0 0. 100-0 0 
080071 6 7 190 4750. 3.2 <13. 0. OO. 0. 0 -1 3 00-0 O- 10 0.0 0 
080081 7 20 00 4750. .538 25. 47. .2373 6.8 0 -1 02 O04. O- 100.0 0 
080091 4 8 0 0 4750. 3.67 11.90 0. O. O~- 1 =-1 3 00-0 O- 100.0 0 
080101 8 9 -2 0 4750. 3.67 -13. 9. 0. O O-l 3 0 O07 O- 190 0.0 0 
080111 9 102 0 4750-5. 3.2 #13. 0. Oj. O- O -1 30056 Of; 100.0 0 
080121 10 19 0 0 4750. ©538 25. 47. .2373 6.8 0 -1 02 0. 0. 1090.0 0 
080131 4190 0 10041.5 3.112 25. 30. 2.2084 1.2.17 2 -1 02 0. O- 19 0. 9 
O80T4l 23 12 0 0 26639. 81-09 17.92 0. 0. O-. 0 =-1 3 00. O- 1091-0 0 
080151 12 13 0 0 26639. 81.09 20.467 0. 0. 9. 0-1 300. O 1901.00 
080161 13 14 0 0 26639. 81.09 232-42 QO. 0. 0. 0 -1 3 0 0. Q. l 0 1.0 0 
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TABLE 6. (cont. ) 
O80171 [4 15 0 O 26639. 81.09 24.42 0. O. 0. O -l 3 OD. O. 1O1-0 D 
080181 15 16 0 O 26639. 81.09 27.17 0. O0O- Oj. O -1 3 OD. O-. LO1.0 OD 
O80191 16 1 0 O 26639. 37-13 29.92 0. O. O. 0 -l1 3 ODO. Of. 101.0 0 
080201 Ll 17 0 0 2944. 2.246 17.0 QO. O- O Ll -1 3 O 1.0 0. 100-0 0 
080211 17 18 0 O 2944. 26.20 18.02 0. Ow. O. O -l1 3 OD. O- 1 ODO 0.0 OD 
080221 18 1 00 2944. 43.635 30.02 0. 0. 0. 0 -1 3 00. O- 1024-0 0 
080231 19 11 00 14791.5 19.69 10.12 0. O. Oj. Ll -l 3 00. 0. 100.0 O 
080241 20 11 0 0 14791.5 19.69 10.12 0. O. O. 1 -1 300. 0 10 7.0 0 
080251 21 5 O11 4750. 3.67 11.0 0. 0. 0. O -1 3 0 0. O 1090 0.0 0 
USOZel 322. 0 30s 484 55-0 le le le O -l OO ODO. O.«. 1 ODO 0.0 O 
080271 21 22 04 O. 3.67 ll. O. le «5 0.22 00. .9 100.0 0 
080281 22 5 04 0. 3.67 ll. O%serea le OD 2 270° 6. 9 10 0-0-0 
080291 11 23 0 O 26639. 20.32 17.0 0. Oj. 90 1-1 300. 0. LOt1-0 0 
080301 0 4 1 O 4074.3135 4.712 40. O. 0. O. O -l1 3 OO. 0. 1900.0 0 
O80311 03 2 2 ~4074.3135 14.75 72. 0. 0. Oj. O -1 3 900. DO 100.0 O 
080321 01300. 1.0 30. O. O. O/. 0 -l 3 OOD. DO. 1 DO O.0 DO 
080331 0 10 4 0 0. 1.00 2 D6 Le 0. 0. 0 -1 0 0 0. 0. l 0 0.0 0 
ROOXXY PUMP NESCRIPTION DATA CARDS 
O90011 3 4 OQ 1 O 1668. 1.0 44943. 710.3 22200. 20000. G7eli 0. 0 
090021 3 4 O 1 O 1668. 1.0 44943. 710.3 22200. 20000. 47.17 O. 0 
O9L001 -11 0. O. ol 0. e1l5 05 024 8 3 - 96 04 - 98 6 97 
091002 8 .9 9 8 96 25 1.0 O. 
092001 -2 0. O. Ll. O. 
O9501L OO. OO. OO. 
095021 OO. 0. 0. 
100000 0 016 0O 
103011 l l 5 0.0 1.31 0.25 1.25 0.5 1.290 0.75 1.12 1.0 1.0 
103021 LL o2.-5°0'«0 -~0.9 025 ~.5 05 -.05 275 4 1.0 1.0 
103031 l 3 5 -1.0 2.1 ~.75 1.775 -0.5 L. 94 -.25 1.4 0. 1.3l 
103041 l 4 5 -1.0 2.1 -.75 1.50 -0.5 1.13 -.25 92 0. 825 
103051 l 5 5 0. 04 025 253 5 625 275 - 80 1.0 1.0 
103061 Lott 5 20% 825 025 ~79 05 815 075 ~90 1.9 1.90 
103071 l 7 #5 -1.0 -1.8 -.75 -l. -.5 -.4 -.25 el 0. 04 
103081 l 8 5 -1.0 -1.8 -.75 -1.72 -.5 -1.54 -.25 -1.27 0. -.9 
L03091 2 l 5 0. 055 025 63 5 71 75 - 86 1.0 1.0 
lLO3S10l 2 2 #5 0. -~.55 02 -.175 4 9 75 -58 1.0 1.0 
LO3LII1 2 3 5 -1.0 1.77 -.75 1.23 -.5 83 -.25 - 62 0. 055 
LOgI2) 2 G@ 5 =1.0 1.77 -.75 1.52 -.5 1.32 -.25 1.13 0. 1.0 
103131 2 5 5 0. -.8 0.2 ~.375 0.5 -.el «75 Pa ie LO ~43 
EOS04) -2° 1°65 70. 1.0 025 88 5 ~75 075 61 1.0 43 
103151 2 7 5 el. -3.45 <=-.75 ~2.75 ~.50 ~-2.0 -.25 -1.375 Q. -.8 
103161 2 8 5 -l. -3.45 -.75 ~2262 --50 -1.583 -.25 -1.15 @Q. -.55 
*LIXXXO VALVE DATA CARDS 
110010 4 0.00. OO. O. OQ. 
110020 5 2 0 0. O. Of OQ. 
110030 -9 OO. O OQ. 0. OO. O. 
110040 +3 0.0900. 0. O. O. 
*L2XXYY LEAK TABLE DATA CARDS 
I2ZOlLOL -2 9 30. Ol. 1000. l. 
120201 6 5 14.7 O«. le 2625 le 2-65 «667 2275 «333 3. O/ LOO00. OQ. 
*LIXXYY FILL TABLE DATA DARDS 
*3XXOO TRIP TYPE POINTS ICALC UNITS PORT HOR X 
1391090 2 I 4 l LBS/SEC 1050. 401.11324 
130101 0. 864.6676 1. 864.6676 4.0 0. 1LOOOO. O. 
1302900 2 2 4 3 LBS/SEC 548. 2999 
130201 OO. Of 940. O. 1025. =-276.2246 3000. -276.2246 
130300 7 2 Ll 3 GAL/MIN 212. 0.0 
130301 0. 15150. 59. 14200. 137. 12500. 178. 11550. 237. 9470. 
130302 268. 7576-4. 280. 5682. 287. 3788. 297. 1894. 304. 0. 1LOOOO. O. 
130400 7 2 9 3 GAL/MIN 212. 0. 
130401 17.6 46355. 23-7 46309. 26.1 4291. G1-2 4175. 42.4 4166. 
130402 98. 3708. 102. 3672. 310. 0. 3000. O. 
*14000 KINETICS CONSTANTS 
*14000 PWR MUL BIL REC UDUF 
140000 3 0 116. 0. l. 
*141001 SCRAM TABLE PATA CARDS 
141001 -8 4 0. 0. 0.9 OO. 1.075 =~.4 1.6 ~4. 2.7 -10. 5.7 -28.1 8. =-30. 10. =-30. 
*142001 DENSITY REACTIVITY DATA CARDS 
142001 -10 .2 -11.55 -3 -8. 830 4 -6.690 -5 -5-Ol °~6 +3.64 e7 -2-51l1 
142002 .8 1.55 °9 +0.73 l. 0. L.1l Q. 
*143001 DOPPLER TA8LE DATA CARDS 
143001 -6 0. 1.85 1000. 9. 2000. -1.71 3000. =-3. 4000. =-4.28 5000. -5.43 





TABLE 6- (cont.-) 


*140XX0 REACTIVITY COEFFICIENT DATA CARDS 

140010 .148289 .148289 OO. O. 

140020 .341782 .341782 0. OQ. 

140030 .019858 .019858 QO. O. 

140040 .341782 .341782 OQ. O. 

140050 .148289 .148289 OQ. QO. 

RYSXXX HEAT SLAB DATA 

ESOC WE SOF 32-00) G0 0 0. 687. 422. “O. 0. 0. 0. OO. Oe 

150012 0. O. 

150021 9 2 20090000. 1368. 698. O. 0. 90. O. O. 0. 

150022 0. QO. 

150031 042 00 0 0 0. 1916. 9772. QO. Of Of DO- DO. Q- 

150032 0. 0. 

50041 0 Tl 2:0 0 0 0 0. 687. 385. 0. 0. O. O. 0. O. 

150042 0. O.- 

150051 4 18 3 000 0 1355. 1355. 239. 0. O. O. DO. O- O. 

150052 0. O. 

150061 124000 0 2288. 2288. 286. 0. O. O. O. O. O.~ 

150062 O. 10.42 

Ls00 71 11 23 5 0000 197. 197. 65. O. 0. 0. 0. O- O. 

150072 0. O. 

150081 16 150000 94. 94. 3le OO. O. 0. GO. O. OW 

150082 0. O. 

150091 11 17 6 000 0 6884. 6884. 95. O Of O- Oe Of O- 

150092 O. 17-61 

150101 0 11 70000 O. 736. 46. O-. O. O/. O- OO O- 

150102 O. O. 

150111 2 3 8000 0 4380. 4380. 184. QO. 0. 0.,0-. O- O- 

150112 0. O. 

50 b2 0 5 905.0 -0070 450. 36. OO. Of O- O. DO. O- 

150122 0. O. 

15053108 9 0 0-0 0 0. 450. 36. (0. 90. 0. 0. 0. 0. 

150132 0. Q-. 

IS014) 079 00 00 0. 838. 67. 0. 0. 0. 0. 0. 0. 

150142 0. O- 

15015) 0 10 9 0000 0, 838. 67. 0. O. Of O- ODO. OD. 

150152 0. QO. 

150 roi 127518 10°60 6 0 1096.8 1096.8 25.67 0. 0. 0. 0. 1Le375 1.375 
150162 O. 2-74 

150171 13 18 10 100 0 1096.8 1096.8 25.67 0. OO. Of Of 446125 4.2125 
PS0172 2.76 5.49 

150181 14 18 10 100 0 398.84 398.84 9.334 0. O. O- OO. 6-0 6-0 
150182 5.51 6.49 

PaO Por t3 16 10 19 0 0 1096.8 1096.8 25.67 9. 0. 06 0. 76875 7.875 
150192 6.51 9.24 

150201 16 18 10 1 0 0 0 1096.8 1096.8 25.67 0. 0. GO. 0. 10.625 10.625 
150202 9.26 11.99 

15021) 0 12 100 0 0 0. 15175.52 178.039 0. .048 0. .0588 1.5375 1.375 
15022 O« 0s 

S022 5-0 13 2.10 0 6 0. 151 752521782039 6. .048 0.00588 45125. 4.125 
Is0222 0. 0. 

150231 01411000 Q0. 5518.372 64.741 90 .048 0. .0588 6.0 6.0 
150232 0. OQ. 

502410 15 1 100 0 0. 5175.52 178.039 0. .048 0. .0588 7.875 7.875 
150242 QO. OQ. 

0251° 0.16.1 1:0 0°00 0. 15175.52) 178.039 0. «048 O. .0588 10.625 10.625 
i>0252 0. 0% i 


*23456789012 3456789012 34567890123456789012 345678901234567890 
*16000 CORE SECTION DATA 


*16000 SLB NODE CLAD OFRAC 
FegO10 21° 7 9 15 0. 017721 
Ieoo20° 22° 7 9 15 0. 26902 
160030523 7 9 15 0. ©10754 
1600490 24 79 15 0. «26902 
Peg0s0725 7 9 15 O- 017721 





*L7TXXYY 


170101 
170102 
170103 
170201 
170202 
170301 
17040} 
170501 
170601 
17070} 
170801 
170901 
171001 
180100 
180101 
180102 
180103 
180200 
180300 
180301 
180302 
180303 
180304 
180305 
180401 
180502 
190100 
190101 
190102 
190103 
190104 
190105 
190200 
190300 
190301 
190302 
190401 
190501 
190502 
200101 
200201 
200301 
200401 
200501 


BR] 


TABLE 6. (cont. ) 
SLAB GEOMETRY DATA CARDS 
2 3 #1 60.0 -920071 1.0 
1 Rv? -000721 0. 
0 3 8 -00267 0. 
1 2 4 3 0.0 04 0%. 
0 5 1 -0208.». 0. 
Wt Ss i120. Asie, GO: 
We ec aes Gee) P25 a, 
Delon? 3 oO 233 O 
ne Gees ee Gan) 20158 0. 
Lee ).eS ele -0625 O. 
l PS BOs 042 0. 
l Pes tO. 08 0. 
1 [x3 le 0234 0. 

-20 500. 3.341 650. 2.971 800. 2-677 950. 724995 TiWO0S 22242 
1250. 2.078 1400. 1.940° 1550. 16823 1700. 1.724 1850. 1.639 
2000. 1.568 2150. 1.507 2300. 1.6457 2450. 1-415 2600. 1.382 
3100. 1.323 3600. 124333 4100. 1.406 4600. 1.538 5100. 1.730 

2 32. 41562 5400. ~41562 
18 32. 7eel2 2 2r6 7.992 392. 8.208 
S72: 8.784 752° 9.540 S32. 10.404 
PEt2s 11.268 [292 . 12.492 1472. 13.176 
1652. 13.968 1832.6 14.796 ZO 2 16.128 
2192, 17.784 2277. 19.656 2552. 21.780 
27,325 24.048 3092. 28.908 3360. 90-120 
HS) Cols O0ca212s 29.5 992s. 28.0 05726-2626 7524 24.7 
-2 200. 8.33 1200. 12.92 
16 B22 3445 122. 38.35 PAP 40095 
392. 43.55 752. 46.8 2012 2-51 35 
C12 we 2 «09 3092. 56655 . 3452. 63.05 
3812. 7268 4352. $9.7 6532. 94.25 
4712. 98.15 4892. !00.1 5144. 101.4 
8000. 101.4 
Z 325 ~-900075 5400. -000075 
5 0.0 28.392 
1480.3 34.476 1675.0 85.176 
L7a755 34.476 3500.0 34.476 
-7 130. 56.9 350. 60-8 450. 62.3 530. 65.2 620. 67.2 710. 70.2 800. 77.5 
-10 68. 52.8 200. 56.7 400. 61-6 600. 64. 800. 66. 1000. .67 
1200. 68.4 1400. 71.8 1600. 75.8 1800. 80.6 
-2 0. 3-718E-6 5000. © 2653E-5 
=2.0.. 0¥ 50004. 0. 
-4 0. 3.0945-6 1652. 4.706E~6 1653. 5.389E-6 5000. 5.389E~6 
=2- "0... 0. 5000. 0. 
=2 0. 0. S000. 0. 
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APPENDIX B 


PLOT/RESTART DATA RETRIEVAL PROGRAM 


This appendix describes a Fortran IV program written to 
extract data from the RELAP4 Plot/Restart tape. During 
RELAP4 problem execution there are two types of information 
written to the Plot/Restart tape. They are called common 
block records, and plot records. Common block’ records 
contain information required to restart the problem. Plot 
records contain data that may be useful for graphical 
display of the transient parameters of interest. 

A plot record is written to the tape at time intervals 
specified by the user (ie. each time aminor edit is 
printed, a plot record is written to the Plot/Restart tape). 
The length of a plot record depends of the system heingp 


modeled and can be calculated as follows, 


LEN = 21 + 20 + 24(NVOL) + 16(NJUN) + 20(NSLB) 


where, 


LEN = length of plot record in computer words, 


NVOL = number of volumes in the RELAP4 model, 


NJUN = number of junctions in the RELAP4 model, 


NSLB = number of heat slabs in the RELAP4 model. 
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For the problem used in this study (ie- Hope Creek) 
NVOL = 23, NJUN = 33, NSLB = 25, and the the plot record 
length (LEN) is thus 1621 computer words. A computer word 
is four bytes or 32 alpha-numeric characters. Future 
reference to word will imply computer word, unless otherwise 
indicated. 

The following is an explanation of the reason for the 
numbers used in the above equation; The first word of each 
plot record contains the letters ‘PLOT’. The next twenty 
words of the plot record contain the problem title and the 
number of junctions and volumes in the system. The next 
twenty words are twenty System parameters, such as 
normalized power, total energy stored in the fuel, transient 
time, etc. (These parameters are listed on pape TI=-228 of 
the RELAP4 Users Manual). The next 24 X NVOL words (in this 
case 552 words) contain the 24 volume parameters listed on 
page II-229 of the Users Manual. This information is 
written to the tape as follows. Words 1-23 are the average 
pressures in volumes 1-23, words 24-46 are the total mass in 
volumes 1-23, etc. The next 16 X NJUN words (in this case 
528 words) contain the 16 junction parameters listed on pape 
II-230 of the Users Manual. This information is written in 
the same order as the volume information (fe- words 1-33 are 
the mass flow rates for junctions 1-33, etc). The next 20 X 
NSLB words (500 words in this case) contain the twenty heat 


slab parameters listed on pages [1-230 and II-23!1 of the 
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Users Manual. These parameters are written in the same 
order as the volume and junction data. 

Table 7 is a listing of the program written to extract 
plotting data from the Plot/Restart tape. The numbers tin 
the left column are sequence numbers used in the following 
program description- The program utilizes a Penn State 
library subroutine called PRECL/FBUF (21). 

The name of the Plot/Restart tape from which data is to 
be extracted must be entered on lines 55 and 1250 (in the 
example shown the tape name is EDIT2). The data set name 
(DSN) on line 1250 must be the same as the as the PSN used 
to write information to the tape during problem execution 
(in the example DSN=TPOUT). 

The first step of the program (line 400) ts the CALL 
PRECL statement. This statement loads the entire first 
record from the Plot/Restart tape into a logical array 
called REC. The length of the record is also’ returned as 
the variable called LEN. A DO loop is then entered (line 
450) where the following steps are taken. Line 500 checks 
the first four characters of the array REC. If the 
characters are not the letters ‘PLOT’ then the record is not 
a plot record, and execution skips to line 650, which is the 
CALL FBUF statement. This statement loads the next record 
of the tape into the array REC and execution returns to line 
450. This process is continued until a plot record is found 


by line 500. 





00050 
00055 
00100 
00150 
00200 
00250 
00300 
00350 
00400 
00450 
00500 
00550 
00600 
00650 
00700 
00750 
00800 
00850 
00900 
00950 
01000 
01005 
01010 
01050 
01100 
01150 
01200 
91250 
01300 
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TABLE 7. PLOT/RESTART TAPE DATA RETRIEVAL PROGRAM 


/x*USERID MVCO1 
/*TAPE EDIT2 
// EXEC FGCG 
//SYSIN DD * 
LOGICAL*1 CNE 
INTEGER*2 LEN 
LOGICAL*1 REC(32000) 
NREC=1 
CALL PRECL(’IN * LEN, NREC,REC) 
DO 1 1#1,9999999 
IF(CNE(REC, PLOT’ ,4))GO TO 1 
KLEN#=LEN /4 
CALL OUT (REC, KLEN) 
1 CALL FBUF(&99,&88) 
99 STOP 
88 STOP 
END 
SUBROUTINE OUT(REC, KLEN) 
REAL REC(KLEN) 
WRITE(17,10)REC(39),REC(52),REC(620),REC(621),REC(622), 
SREC (229) ,REC(1169),REC(1219) 
WRITE(6, 10) REC(39),REC(52),REC(620),REC(621),REC(622), 
SREC(229),REC(1169),REC(1219) 
10 FORMAT (F8.3, 1P7E10-3) 
RETURN 
END 
//DATA-FT17F001 DD UNIT#=BAT, FILES=$DEM6 
//DATA-IN DD UNIT=2400, VOL=SER#=EDIT2, DSN=TPOUT 
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When a plot record is found program execution goes to 
line 550 which divides the record length LEN by four, 
converting the record length from computer words’ to alpha- 
numeric words. The subroutine OUT is then called (line 600) 
and execution skips to line 850. Line 900 changes the array 
REC from a logical array to a real number array. Lines 959 
and 1000 write the desired parameters (discussed helow) to 
Fortran Unit 17, defined in line 1200 as a BAT file called, 
in this case SDEM6. Lines 1005 and 1OLO write the same 
parameters to Fortran Unit 6, a line printer. Following 
printing of the desired parameters, program execution is 
returned to the main program at line 650. The program 
continues to process”7 records until an end-of-file is 
encountered, when execution is stopped. 

To determine the word number of the desired parameter 
to be printed, the following steps are taken; 

For system parameters, the word number is simply 21 
plus the number of the systen parameter. For example, the 
transient time is parameter number 18 (recall that parameter 
numbers are taken from the RELAP4 Users Manual on the papes 
noted earlier). Thus the word number for transient time is 
21 + 18 = 39 and REC(39) contains this parameter. 

For volume parameters, the word number is calculated as 


follows, 


WN = 41 + (VPN - 1)NVOL + VOLN 
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where the new symbols are, 


WN = word number of desired parameter 

VPN = volume parameter number 

VOLN = volume number of desired parameter. 

For example, if the desired parameter is mixture level 
in volume 4, tthe word number is calculated as follows; The 
VPN for mixture level is 9. NVOL = 23, and VOLN = 4. Thus 
WN = 229, and REC(229) is the mixture level in volume 4. 

For junction parameters the word number tis calculated 


from, 


WN = 41 + 24(NVOL) + (JPN - L)NJUN + JUNN 


where the new symbols are, 


JPN = junction parameter number 

JUNN = junction number of desired parameter. 

For example, if the desired parameter is the flow rate 
in junction 27, the word number is calculated as follows; 
The JPN for flow rate is !, NVOL = 23, NJUN = 33, and JUNN = 
27. Thus WN = 620 and REC(620) is the flow rate at junction 
27. 

For heat slab parameters the word number is calculated 


from, 


WN = 41 + 24(NVOL) +16(NJUN) + (SLBPN = 1)NSLB + SLBN 








RR 


where, 
SLBPN = heat slab parameter number 
SLBN = heat slab number. 

For example, if the desired parameter is surface 
temperature, heat slab $23, the word number is calculated as 
follows; The SLBPN for surface temperature is 2, NVOL = 23, 
NJUN = 33, NSLB = 25, and SLBN «= 23. Thus WN = 1169, and 
REC(1169) is the surface temperature of heat Slab 23. 

The central processing unit (CPU) time required for 
execution of the program depends on the amount of data to be 
processed- As a benchmark, processing of 400 records takes 


about fifty CPU seconds. 
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